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INTRODUCTION. 


$1. A great many if not all of the problems in mathematics 
may be so formulated that they consist in finding from given 
data the values of certain unknown quantities subject to certain 
conditions. We may distinguish different stages in the solution 
of a problem. The first stage we might say is the proof that the 
quantities sought for really exist, that it is possible to satisfy 
the given conditions or, as the case may be, the proof that it is 
impossible. In the latter case we have done with the problem. 
Take for instance the celebrated question of the squaring of the 
circle. We may in a more generalized form state it thus: Find 
the integral numbers, which are the coefficients of an algebraic 
equation, of which 7 is one of the roots. Thirty years ago 
Lindemann showed that integral numbers subject to these con- 
ditions do not exist and thus a problem as old almost as 
human history came to an end. Or to give another instance 
take Fermat’s problem, for the solution of which the late Mr. 
Wolfskehl, of Darmstadt, has left $25,000 in his will. Find the 
integral numbers 2, y, z that satisfy the equation 


am + De —< ie 


where 7 is an integral number greater than two. Fermat main- 
tained that it is impossible to satisfy these conditions and he is 
probably right. But as yet it has not been shown. So the 
solution of the problem may or may not end in its first stage. 

In many other cases the first stage of the solution may be so 
easy, that we immediately pass on to the second stage of finding 
methods to calculate the unknown quantities sought for. Or 
even if the first stage of the solution is not so easy, it may be 
expedient to pass on to the second stage. For if we succeed in 
finding methods of calculation that determine the unknown quan- 
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tities, the proof of their existence is included. If on the other 
hand, we do not succeed, then it will be time enough to return 
to the first stage. 

There are not a small number of men who believe the task of 
the mathematician to end here. This, I think, is due to the 
fact that the pure mathematician as a rule is not in the habit of 
pushing his investigation so far as to find something out about the 
real things of this world. He leaves that to the astronomer, to 
the physicist, to the engineer. These men, on the other hand, 
take the greatest interest in the actual numerical values that 
are the outcome of the mathematical methods of calculation. 
They have to carry out the calculation and as soon as they do so, 
the question arises whether they could not get at the same result 
in a shorter way, with less trouble. Suppose the mathematician 
gives them a method of calculation perfectly logical and con- 
clusive but taking 200 years of incessant numerical work to 
complete. They would be justified in thinking that this is not 
much better than no method at all. So there arises a third stage 
of the solution of a mathematical problem in which the object is 
to develop methods for finding the result with as little trouble as 
possible. I maintain that this third stage is just as much a 
chapter of mathematics as the first two stages and it will not do 
to leave it to the astronomer, to the physicist, to the engineer or 
whoever applies mathematical methods, for this reason that 
these men are bent on the results and therefore they will be apt 
to overlook the full generality of the methods they happen to 
hit on, while in the hands of the mathematician the methods 
would be developed from a higher standpoint and their bearing 
on other problems in other scientific inquiries would be more 
likely to receive the proper attention. 

The state of affairs today is such that in a number of cases the 
methods of the engineer or the surveyor are not known to the 
astronomer or the physicist, or vice versa, although their prob- 
lems may be mathematically almost identical. It is particularly 
so with graphical methods, that have been invented for definite 
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problems. A more general exposition makes them applicable 
to a vast number of cases that were originally not thought of. 

In this course I shall review the graphical methods from a 
general standpoint, that is, I shall try to formulate and to teach 
them in their most generalized form so as to facilitate their 
application in any problem, with which they are mathematically 
connected. The student is advised to do practical exercises. 
Nothing but the repeated application of the methods will give 
him the whole grasp of the subject. For it is not sufficient to 
understand the underlying ideas, it is also necessary to acquire a 
certain facility in applying them. You might as well try to learn 
piano playing only by attending concerts as to learn the 
graphical methods only through lectures. 

1 For the literature of the subject see “Encyklopidie der mathematischen 


Wissenschaften,” Art. R. Mehmke, ‘Numerisches Rechnen,” and Art. F 
Willers and C. Runge, ‘‘Graphische Integration.” 
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CHAPTER I. 
GRAPHICAL CALCULATION. 


§ 2. Graphical Arithmetic.—Any quantity susceptible of mensu- 
ration can be graphically represented by a straight line, the 
length of the line corresponding to the value of the quantity. 
But this is by no means the only possible way. A quantity 
might also be and is sometimes graphically represented by an 
angle or by the length of a curved line or by the area of a square 
or triangle or any other figure or by the anharmonic ratio of four 
points in a straight line or in a variety of other ways. The 
representation by straight lines has some advantages over the 
others, mainly on account of the facility with which the ele- 
mentary mathematical operations can be carried out. 

What is the use of representing quantities on paper? It is a. 
convenient way of placing them before our eye, of comparing 
them, of handling them. If pencil and paper were not as cheap 
as they are, or if to draw a line were a long and tedious under- 
taking, or if our eye were not as skillful and expert an assistant, 
graphical methods would lose much of their significance. Or, 
on the other hand, if electric currents or any other measurable 
quantities were as cheaply and conveniently produced in any 
desired degree and added, subtracted, multiplied and divided 
with equal facility, it might be profitable to use them for the 
representation of any other measurable quantities, not so easily 
produced or handled. 

The addition of two positive quantities represented by straight 
lines of given length is effected by laying them off in the same 
direction, one behind the other. The direction gives each line a 
beginning and an end. The beginning of the second line has to 
coincide with the end of the first, and the resulting line represent- 
ing the sum of the two runs from the beginning of the first to 
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the end of the second. Similarly the subtraction of one positive 
quantity from another is effected by giving the lines opposite direc- 
tions and letting the beginning of the line that is to be subtracted 
coincide with the end of the other. The result of the subtrac- 
tion is represented by the line that runs from the beginning of 
the minuend to the end of the subtrahend. The result is positive 
when this direction coincides with that of the minuend, and nega- 
tive when it coincides with that of the subtrahend. This leads 
to the representation of positive and negative quantities by lines of 
opposite direction. The subtraction of one positive quantity from 
another may then be looked upon as the addition of a positive and 
a negative quantity. I do not want to dwell on the logical explana- 
tion of this subject, but I want to point out the practical method 
used for adding a large number of positive and negative quantities 
represented by straight lines of opposite direction. Take a 
straight edge, say a piece of paper folded over so as to form a 
straight edge, mark a point on it, and assign one of the two 
directions as the positive one. Lay the edge in succession over 
the different lines and run a pointer along it through an amount 
equal in each case to the length of the line and in the positive 
or negative direction according to the sign of the quantity. The 
pointer is to begin at the point marked. The line running from 
this point to where the pointer stops represents the sum of the 
given quantities. The advan- 
tage of this method is that the 
intermediate positions of the 
pointer need not be marked pro- 
vided only that the pointer keeps 
its position during the move- 
ment of the edge from one line 
Fie 1. to the next. As an example take 

the area, Fig.1. A number of 

rectangular strips } cm. wide are substituted for the area so that, 
measured in square centimeters, it is equal to half the sum of 
the lengths of the strips measured in centimeters. The straight 
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edge is placed over the strips in succession and the pointer is 
run along them. The edge is supposed to carry a centime- 
ter scale and the pointer is to begin at zero. The final position 
of the pointer gives half the value of the area in square centi- 
meters. ‘The drawing of the strips may be dispensed with, their 
lengths being estimated, only their width must be shown. If 
the scale should be too short for the whole length, the only thing 
we have to do is to break any of the lengths that range over the 
end of the scale and to count how many times we have gone 
over the whole scale. I have found it convenient to use a little 
pointer of paper fastened on the runner of a slide rule so that it 
can be moved up and down the metrical scale on one side of the 


slide rule. The area is in this manner determined rapidly and 
with considerable accuracy, very well comparable to the ac- 
curacy of a good planimeter. If the area of any closed curve 
is to be found, the way to proceed is to choose two parallel 
lines that cut off two segments on either side (see Fig. 2), to 
measure the area between them by the method described above 
and to estimate the two segments separately. If the curves of 
the segments may with sufficient accuracy be regarded as arcs 
of parabolas the area would be two thirds the product of length 
and width. If not they would have to be estimated by substitut- 
ing a rectangle or a number of rectangles for them. 
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In the same way the addition and subtraction of pure numbers 
may also be carried out. We need only represent the numbers 
by the ratios of the lengths of straight lines to a certain fixed 
line. The ratio of the length of the sum of the lines to the length 
of the fixed lines is equal to the sum of the numbers. The con- 
struction also applies to positive and negative numbers, if we 
represent them by the ratio of the length of straight lines of 
opposite directions to the length of a fixed line. 

In order to multiply a given quantity ¢ by a given number, 
let the number be given as the ratio of the lengths of two straight 
lines a/b. If the quantity c is also represented by a straight line, 
all we have to do is to find a straight line 2 whose length is to 
the length of c as a to b. This can be done in many ways by 


---— 


Fra. 3. Fig. 4. 


constructing any triangle with two sides equal to a and b and 
drawing a similar triangle with the side that corresponds to 6b made 


equal to c. Asa rule it is convenient to draw a and b at right \ 


angles and the similar triangle either with-its hypotenuse parallel 
(Fig. 3) or at right angles (Fig. 4) to the hypotenuse of the first 
triangle. Division by a given number is effected by the same con- 
struction; for the multiplication by the ratio a/b is equivalent 
to the divisions by the ratio b/a. 

If a, b, ¢ are any given numbers, we can represent them by the 
ratios of three straight lines to a fixed line. Then the ratio of 


ef 
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the line constructed in the way shown in Fig. 3 and Fig. 4 to 
the fixed line is equal to-the number 


ac 
b* 

Multiplication and division are in this way carried out simul- 
taneously. In order to have multiplication alone, we need only 
make 6 equal | and in order to have division alone, we need only 
make a or ¢ equal 1. 

In order to include the multiplication and division of positive 
and negative numbers we can proceed in the following way. Let 
the lines corresponding to a, 2, Fig. 3, be drawn to the right side 
of the vertex to signify positive numbers and to the left side to 
signify negative numbers. Similarly let the lines corresponding 
to b, c be drawn upward to signify positive numbers and down- 
ward to signify negative numbers. Then the drawing of a 
parallel to the hypotenuse of the rectangular triangle a, b through 
the end of the line corresponding to e will always lead to the 


number 


ac 
t= 
b 


whatever the signs of a, b, c may be. 

The same definition will not hold for the construction of Fig. 4. 
If the positive direction of the line corresponding to a is to the 
right and the positive direction of the line corresponding to 6 is 
upwards then the positive directions of x and ¢ ought to be such 
that when the right-angled triangle x, ¢ is turned through an 
angle of 90° to make the positive direction of zx coincident 
with the positive direction of a, the positive direction of ¢ coin- 
cides with the positive direction of 6. If we wish to have the 


positive direction of x upward, the positive direction of ¢ would 


have to be to the left, or if we wish to have the positive direction 
of c to the right, the positive direction of x would have to be 
downward. If this is adhered to, the construction for division 
and multiplication will include the signs. 
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§ 3. Integral Functions.—We have shown how to add, subtract, 
multiply, divide given numbers graphically by representing them 
as ratios of the lengths of straight lines to the length of a fixed 
line and finding the result of the operation as the ratio of the 
length of a certain line to the same fixed line. By repeating 
these constructions we are now enabled to find the value of any 
algebraical expression built up by these four operations in any 
succession and repetition. Let us see for instance how the values 
of an integral function of 2, that is to say, an expression of the form 


Go + ax + apa? + --- + aya” 


may be found by geometrical construction, where do, a +++ Gn, & 
are any positive or negative 
numbers. We shall first as- 
sume that all the numbers are 
positive, but there is not the 
least difficulty in extending 
the method to the more gen- 
eral case. 

Now let ao, a1, 2, +++ An 
signify straight lines laid off 
on a vertical line that we call 
the y-axis, one after the other 
as if to find the straight line 


0 x 1 oy 


Fic. 5. io + di -F Ga > a 


The lengths of these lines measured in a conveniently chosen 
unit of length are equal to the numbers designated by the same 
letters. In Fig. 5 ao runs from the point O to point Ci, a from 
Cy, to Cy, --- a, from C, to Crs. 

Let x be the ratio of the lines Ox and 01, Fig. 5, drawn hori- 
zontally from O to the right. The length O1 is chosen of con- 
venient size independent of the unit of length that measures the 
lines do, a1, ++ Gn. The length Oz is then defined by the value 
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of the ratio a. Through a and 1 draw lines parallel to the y-axis. 
Through C,,4: draw a line parallel to Oz, that intersects those two 


parallels in P,, and B,. Draw the line B,,C,, that intersects the 


parallel through 2 in P,4.. Then the height of P,; above C, 
will be equal to a,x. For if we draw a line through P,,_; parallel 
to Ox intersecting the y-axis in D,, the triangle C,D,P,-1 will be 
similar to C,zCniiB, and their ratio is equal to 2, therefore 
C,D, = anx. Consequently the height of P,; above Cy: is 
equal to Cyr1D, = anz + dn. Now let us repeat the same 
operation in letting the point D, take the part of Cry. Through 
D,, draw a line parallel to Oz, that intersects the parallels through 
xand lin Py; and B,_;. Draw the line B,_1C,_1 that intersects 
the parallel through zx in P,_». 
Then the height of Pn». above 
C,-1 will be equal to 


Cy aD,-2 = (G,2 + Gy-1)2, 


and the height above C,_2 will be 
equal to 


ar + An—1 at An-2- 


Continue in the same way. Draw 
P,2Bn2 parallel to Oz, draw 
B,-2Cn-2 and find the point P,-3. 
Then the height of P,-3 above Cy_2 will be 

(Gnd + One + On—2) 
and the height of P,-3 above Cn-s 


Ont? + Ant" + Ono + On-3- 


Finally a point Pp is found (see Fig. 6 for n = 4) by the inter- 
section of B,C; with the parallel to the y-axis through x, whose 
height above O is equal to 

Galt te Oy at" of 4 et aie t- Oo 
Let us designate the line ~Po by y, so that 


8 GRAPHICAL METHODS. 


Y = Ont” + dnae™* + +++ + Ue + ao, 


in the sense that y is a vertical line of the same direction and 
length as the sum of the vertical lines an#”, Gna”, +++ GX, do. 
The same construction holds good for values of z greater than 
1 or negative. The only difference is that the point 2 is beyond 
the interval 01 to the right of 1 or to the left of O. The negative 
sign of 
Onl, Ont tly, Ont ees Sle, 


will signify that the direction of the lines is downward. Nor are 
any alterations necessary in order to include the case that several 
or all of the lines ao, a, --+ dn are directed downward and corre- 
spond to negative numbers. They are laid off on the y-axis in 
the same way as if to find the sum 


dye as ES ok ee 


C.41 lying above or below C, according to a, being directed 
upward or downward. The construction can be repeated for a 
number of values of 2. The points Po will then represent the 
curve, whose equation is 


Y = Got Ge +++ + ane", 


x and y measuring abscissa and ordinates in independent units 
of length. 

In order to draw the curve for large values of x a modification 
must be introduced. It will not do to choose 01 small in order 
to keep 2 on your drawing board; for then the lines B,C, will 
become too short and thus their direction. will be badly defined. 
The way to proceed is to change the variable. Write for instance 
X = 2/10, so that X is ten times as small as 2 and write 


A. = ole 
Then as 


y = a+ a-1075+ a2102 


Dod Ay 
102+ -+- + @,10 10" 


he 
> 
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we find 
y = Ao + AVX + AgX? + ++» + ApX 


Lay off the lines Ao, Ai, --- A, in a convenient scale and let 
X play the part that x played before. The curve differs in scale 
from the first curve and the reduction of scale may be different 
for abscissas and ordinates but may if we choose. be made the 
same so that it is geometrically similar to the first curve reduced 
to one tenth. It is evident that any other reduction can be 
effected in the same manner. By increasing the ratio 2/X we 
enhance the value of A, in comparison to the coefficients of lesser 
index, so that for the figure of the curve drawn in a very small 
scale all the terms will be insignificant except A,X”. In this 
case the points Ci, C2, ---, Cy will very nearly coincide with O 
and only Cy41 will stand out. 

It is interesting to observe that the best way of calculating an 
integral function 


Go aye + aa" + --* + a2” 


for any value of x proceeds on exactly the same lines as the 


geometrical construction. The coefficient a, is first multiplied 


with 2 and a,_1 is added Call the result a,_1’._ This is again 
multiplied by x and a,-2 is added. Call this result a,»’.. Con- 
tinuing in this way we finally obtain a value of ao’, which is equal 
to the value of the integral function for the value of « considered. 
Using a slide rule all the multiplications with 2 can be effected 
with a single setting of the instrument. The coefficients a, and 
the values a,’ are best written in rows in this way 


Qn On-1 OUn2 °**° 1 ao 
Gf Unt <> t  Oy'2 

7 

Grr Gyno’ +++ ay’ ao 


The accuracy of the slide rule is very nearly the same as the 
accuracy of a good drawing. But the rapidity is very much 
greater. When therefore only a few values of the integral func- 


tion are required, the geometrical construction will not repay 
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the trouble. It is different, however, when the object is to make 
a drawing of the curve. The values supplied by calculation 
would have to be plotted, while the geometrical construction 
furnishes the points of the curve right away and in this manner 
gains on the numerical method. 

There is another geometrical method, which in some cases 
may be just as good. Let us propose to find the value of an 
integral function of the fourth degree. 


Y = do + aye + ax? + aga? + agt* 


and let all coefficients in the first instance be positive. 
The coefficients ao, a1, 2, 43, dg are supposed to be represented 
by straight lines, while x will be the ratio of two lines. The lines 
lp, Q1, Az, a3, a4 are laid off ina 
broken line ap to the right from 
Co to Ci, a, upward from C; to 
C2, az to the left from C; to C3, a3 
downward from C3 to C4, a4 again 
to the right from C, to C; (Fig. 7). 
Through C; draw a line C;A to 
a point A on C;3C;, or its prolonga- 
‘tion and let x be equal to the 
ratio Cyd : CsCs taken positive 
when (4A has the same direc- 


tion as C;C,. Then we have 


CA = ag, 


and 
C3A = Ase + a3. 


C,A and C3A are positive or negative according to their direction, 
being the same as the direction of C3C, or opposite to it. Through 
A draw the line AB forming a right angle with 0;A to a point B 
on C2C3 or its prolongation. Then we have 


C3B = C3 A-x = (ayx + a3) x 
and 
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CB = ag? + age + aw. 


C3B and CB are positive or negative according to their direction 
being the same as the direction of C.C3 or opposite to it. Simi- 
larly we get 


CLD = aga® + asx? + ae + am, 
and finally 


CoE = agat + aga? + aca? + aye + ao. 


CoE is positive, when F is on the right side of Cy and negative 
when on the left side. When the point A moves along the line 
C3Cs, the point E will move 
along the line C,C; and its 
position will determine the 
values of the integral function. 
To find the position of E for 
any position of A, we might 
use transparent squared paper, 
that we pin onto the drawing 
at C;, so that it can freely be 
turned round C;. Following 
the lines of the squared paper 
along C;ABDE after turning it through a small angle furnishes 
the position of E for a new position of A (Fig. 8). 

To include the case of negative coefficients we draw the corre- 
sponding line in the opposite direction. If for instance ag is 
negative C3C, would have to lie above C3; but C3A would have 
to be counted in the same way as before, positive in a downward, 
negative in an upward direction. 

The extension of the method to integral functions of any degree 
is obvious and need not be insisted on. It may be applied with 
advantage to find the real roots of an equation of any degree. 
For this purpose the broken line C;ABDE would have to be 
drawn in such a way that E coincides with Co. In the case of 
Fig. 7, for instance, it is easily seen that no real root exists. 
Fig. 9 shows the application to the quadratic equation. A circle 


Fie. 8. 
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is drawn over C,(3 as diameter. Its intersections with CC, 
furnish the points A and A’ that correspond to the two roots. 
Both roots are negative in this case. 
The first method of constructing 
the values of an integral function can 
4 be extended to the case where the 
ad be function is given as the sum of a 
number of polynomials of the form 


Q 
ee 


ie 


Y = A+ a(t — p) + a(x—p)(e —g) 
+ ag(z — p)(e — Qe — 1) 4s 


Ch Let us again suppose do, a, de, «°- 
to represent straight lines laid off as 
before on the y-axis upwards or down- 
wards as if to find their sum. 2, p,q, r --- are meant, to be num- 
bers represented by the ratio of certain segments on the axis of 
abscissas. Let us consider the case of four terms, the highest poly- 
nomial being of the third degree. The fixed distance between the 
points marked p and p+1,qand ¢g+ 1, r and r+ 1 on the 
axis of abscissas, Fig. 10 is chosen arbitrarily and the position 


~ 


~ 
~ 


Fig. 10. 


of the points marked p, q, r, x is made such that the ratio of 
Op, Og, Or, Ox to that fixed distance is equal to the numbers 


p, q, 7, *. For negative values the points are taken on the left 
of O. 
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j en 

3 

Draw parallels to the y-axis through p, q, 7, x, p+ 1,q+1, 
r+ 1. On the parallel through r+ 1 find the point Qo of the 
same ordinate as C, and on the parallel through r find the point 
Ao of the same ordinate as C3. Join Ap and Qo by a straight 
line and find its intersection P; or that of its prolongation with 
the parallel through x. The height of P, above C3 or Apo is 
equal to a3(z — r) and the height above C, is equal to a3(a — r) 
+ ad. On the parallel through g + 1 find a point Q, of the same 
ordinate as P; and on the parallel through qg a point A; of the 
same ordinate as C2. Join A; and Q, by a straight line and find 
its intersection P, or that of its prolongation with the parallel 
througha. The height of P; above C; or A; is equal to 


las(z@ — r) + w)](x — q), 
and the height above C; is equal to 
as(x — r)(@ — q) + (x — gq) + G. 


Finally find a point Q on the parallel through p+ 1 of the 
same ordinate as P, and a point A» on the parallel through p of 
the same ordinate as C;. Join A: and Q: by a straight line and 
find its intersection P3 or that of its prolongation with the par- 
allel through z. The height of P; above Ci or A: will then be 
equal to 

fant — 1) (2 — g) + an(z — g) + alle — p) 


and the ordinate of P3 will be equal to the given integral function 


i= Og(t — 1) (x — g)(e@ ~p) + m(x — g)(e — p) 
+ a(x a D) + ao. 


For large numbers 7p, q, r, we use a similar device as before by 
introducing new numbers P, Q, R, X equal to one tenth, or one 
hundredth or any other fraction of pgrz. For instance 


P=p/10, Q=9/10, R=7/10 X =7/10. 
We then write 
Ao = ao, Ay = 10a, Ay = 100a2, Az — 1000az, 
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and obtain 
y = Ap + A(X — P) + A(X — P(X — Q) 
+ A,(X — P)(X — Q)( t= Be 

The scale for the lines Ao, A1, Az, A3 and y must then be reduced 
conveniently and the values are constructed in the same way as 
before. 

Now let us consider the inverse problem. The values of the 
integral function are given for 


re od! ee Pe 
find the lines ao, a1, a2, a3, So that the value of the integral function 
may be found for any other value of 2 in the way shown above. 

Let us designate the given values of the integral function for 
L= p, 9,7; 8 DY Yps Yar Yrs Ys and the points on the parallels through 
Dp; q, ’, § with these ordinates by P, Q, R, S (see Fig. 12). 

For z = p the integral function 


Y = a + a(@ — p) + a(x — p)(x — q) + a3(x— p)(ex—g)(a—1) 
reduces to ao. Therefore we have y» = do. The point C; is 
found Le drawing a parallel to the axis of abscissas through P 
and taking its intersection with 
the axis of ordinates. 

In order to find C, draw a 
straight line through P and Q 
and find its intersection A with 
the parallel through p + 1 (Fig. 
11). A parallel to the axis of 
abscissas through A_ intersects 
the axis of ordinates in (>. a 

Fief 11: the differences yg—yp and Ya— 
(writing y. for the ordinate ot 
A) are proportional to the differences of the abscissas and con- 
sequently in the ratio (g — p):1. Therefore 


Yo Y 
Yai Ye ge ee 


i fs 


+ 
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In the same way as the point Q on the parallel through q we 
might join any point X ona parallel through x with the point P, 
find the intersection with the parallel through p + 1 and draw a 
parallel to the axis of abscissas. The point of intersection of 


Fie. 12; 


this parallel with the vertical through 2 let us call X’ and its 
ordinate y’. Then we have 


ran 
of — tp = = + ale — 9) + ale — De), 
Let us carry out this construction not only for x = q but also 
for x=randx=s. This leads us to three points Q’, R’, S’ 
on the verticals through q, r, s, whose ordinates are the values 


of the integral functions 
y’ = (do + a) + (x — 9) + as(e — g)(@ — 1). 

In this way we have reduced our problem. Instead of having 
to find an integral function of the third degree from four given 
points P, Q, R, S, we have now only to find an integral function 
of the second degree from three given points Q’, R’, S’.. A second 
reduction is effected in exactly the same manner. Q’ is joined 
with R’ and S’ by straight lines and through their intersection 
with the vertical through g + 1 parallels to the axis of abscissas 
are drawn that intersect the verticals through r and s in the 
points R’’ and 8” respectively. The ordinates of these points 
are the values of the integral function y’’ defined by 
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, 


y’ — Ya 
x =F 


y' — Yq = = m+ a(x— 1), 


for « = rand z = 8, or 
y" =a tatat a(x — 1). 

The horizontal through R” intersects the axis of ordinates in the 
point (3. Finally we find C, by drawing a parallel to the axis 
of abscissas through the intersection of R’’S” or its prolongation 
with the vertical through r + 1. 

Having found the points 0:02C;C, we can now for any value 
of x construct the ordinate 


y = a t+ a(x — p) + a(x — p)(@ — g) 
+ as(z — pyle — 9) ae 


and thus draw the parabola of the third degree passing through 
the four points P, Q, R, S. 

The construction may be somewhat simplified first by making 
p+i1=gq. Our data are the points P, Q, R, S, and we are 
perfectly at liberty to make the vertical through p + 1 coincide 
with the vertical through Q. In this case the point Q’ will 
coincide with Q. The parabola of the second degree through the 
points Q’R’S’ is again independent of the distance between the 
verticals through g and g + 1 and at the same time independent 
of the point P. Therefore we are perfectly at liberty, for the 
construction of any point of this parabola, to make the vertical 
through q+ 1 coincide with the vertical through R even if the 
distance of the verticals through P and Q is different from that 
of the verticals through Q and R. R” will in this case coincide 
with R’. The procedure is shown in Fig. 12. Starting from 
the points P, Q, Rf, S the first step is to find R’, 8’ by, connecting 
R and S with P and drawing horizontals through the inter- 
sections A,, A, with the vertical through g. The next step is to 
find S’’ by connecting Q (identical with Q’) with S’ and drawing 
a horizontal through the intersection with the vertical through r. 
Now the straight line R’’S’”” can be drawn (R” being identical 
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with R’). On the vertical through any point x take the inter- 
section with R’’S” and pass horizontally to the point A,’ on the 
vertical through r. Draw the line Q’ 4,’ and find its intersection 
with the vertical through x. This point is on the parabola 
through Q’R’S’. Pass horizontally to the point A, on the 
vertical through g and draw the line A,P. Its intersection with 
the vertical through 2 is a point on the parabola of the third 
degree through P, Q, R, S. 

The method is evidently applicable to any number of given 
points, the degree of the parabola being one unit less than the 
number of points. 

The methods for the construction of the values of an integral 
function may be applied to find the value of any rational function 


y= R@) "| 


For a rational function can always be reduced to the form of a. 
quotient of two integral functions 


R(x) = gi(x)/go(a). 


Now after having constructed curves whose ordinates give the 
values of g:(x) and go(x) for any abscissa x (Fig. 13), R(x) is found 
in the following manner. 
Through a pointP on the 
axis of abscissas draw a 
parallel to the axis of or- 
dinates. Let G; and G 
be the points whose ordi- 
nates are equal to (2) 
and go(x). Pass horizon- 
tally from G; to Gy’ on the Fic. 18, 

vertical through P and 

from G, to Gy’ on the axis of ordinates. Draw a line through 
P and G,' and produce it as far as A where it intersects the 
horizontal through G,. Then R(x) is equal to the ratio G,’A 
to PO. G,'A may then be set off as ordinate on the vertical 

3 
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through 2 and defines the point M whose ordinate is equal to 
R(z) in length, when OP is chosen as the unit of length. 

§ 4. Linear Functions of Any Number of Variables.—Let us 
consider a linear function of a number of variables 24, a «++ 2n, 


Qo + arty + dna, + +++ + An®n, 


where do, a, d2, *** Gn are given numbers positive or negative. 
The question is how the value of this linear function may be 
conveniently constructed for various systems 21, %, +++ 2n- 
Suppose do, a, --- @n to represent horizontal lines directed to 
the right or left according to the sign of the corresponding number 
and to be laid off on an horizontal axis in succession as if to find 
the sum 


Go “+ & + G2 FF =<" - Ga, 


ao begins at O and runs to (Cj, a begins at C; and runs to C, and 
so on (Fig. 14). The numbers 21, a, --- 2, let us represent 


Fie. 14. 


by ratios of lengths. We draw a vertical line through O and 
choose a point P on the horizontal axis. Then let a, be equal 
to the ratio 01/PO, 2 = 02/PO, etc. If P is chosen on the left 
of O, we take the point 1 above O for a positive value of a, and 
below O for a negative one and the same for the other points. 
Mark a point 0 above O in the same distance from 0 as P. Join 
the point P with the points 0, 1, 2, 3, 4, --- and draw a broken 


cia 


of 
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line 0A A142A3Aq in such a manner that Ap is on the vertical 
through C, and Op is parallel to PO, A, on the vertical through 
Cz and AoA; parallel to P1, Az on the vertical through C3 and 4,4, 
parallel to P2 and so on. Then the ordinate yo of Ao will have. 
the same length as ao and will be directed upward when the 
direction of a is to the right, and downward when the direction 
of ao is to the left. The difference y; — yo of the ordinates of A; 
and Ap is equal in length to ma, as y: — yo and a, have the same 
ratio as Ol and PO. A, will be above or below Ap according to 
the line a;2; being directed to the right or to the left and it is 


understood that az; has the same direction as a for positive 


Cn CO oiyAn 
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values of 2; and a direction opposite to a for negative values 
of 2; Thus the ordinate y; has the same length as the line 
ap + a2; and its direction is upward or downward according to 
the direction of the line ap + a2 being to the right or to the left. 
In the same way it is shown that the ordinate y% of the point A>» 
is equal in length to 
Ae + a1 + AX, 
and y3 to 
Qo + aX + Aet, + a3Xs 


and so on, the direction upward or downward corresponding 
to the positive or negative value of the linear function. 
If the values of 2, a, +++ an satisfy the equation 


Gig + aya + det, + +++ + Antn = 0 


the ordinate y, must vanish, that is to say, the point A, must 
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coincide with C41, the end of the line a,. And vice versa if A, 
and Cn41 coincide the equation is satisfied. Consequently if we 
know all the values but one of the numbers 2, 2, +++ a the 
unknown value can be found graphically. For suppose 23 to be 


CriCn Cnei=An 
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the unknown value we can, beginning from QO, find the broken 
line as far as A, and beginning from the other end A, we 
can find it as far as A3 (Fig. 15). A parallel to A:A; through P 
furnishes the point 3 on the axis of ordinates. If a, a, +++ &4 
are known and only 2, not, we can draw the broken line as far 
as An and as A, has to coincide with C,41, we can draw a parallel 
to An+A, through P and find the point n on the axis of ordinates 


O GQ Co Us OsCy C5C7 


Bre. 17. 


that determines the value x, by the ratio On/PO or On/Oo. In 
Figs. 15 and 16 all the coefficients ao, a, ---, are positive. A 
negative coefficient a; is shown in Fig. 17. The only difference 
is that C, lies to the left of C; and consequently the broken line 
passes from A, back to As. 

If we keep the points 0, 1, 2, ---, in their positions but change 
the position of P to P’ (Fig. 18) and repeat the construction of 


i 


i 
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the broken line, we obtain 0A)’ Ay’ A,’ «+ instead of OAjAyA2 +++. 
The ordinate y,’ of the point A,’ is evidently 
Ya Samat ant bape 
and therefore 
Yat 
P’0 


That is to say, by changing the position of P without changing 
the position of the points 0, 1, 2, --- we can change the scale of 
the ordinates of the broken line. They change inversely pro- 


, 


i = 
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portional to PO. It may be convenient to make use of this 
device in order to make the ordinates a convenient size inde- 
pendent of the scale that we have chosen for the coe OQ ie2; 
that determine the values 
ho ears 
Sari 00, 
A linear equation with only one unknown quantity 
ao + aa, = 0 


is solved by drawing a parallel to 4A; through P. Let a second 
equation be given with two unknown quantities 


bo a. bya -- boats =, 


The lines bo, 61, b: are laid off as before. Knowing 2 as the 
solution of the first equation we can construct the broken line 
OB,B1 corresponding to the second equation and as B, must 
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coincide with the end of hk, we can draw a parallel to BiB, 
through P and find a. In a similar manner we can find 23 
from a third equation 


CoH Cyt + Cote + C3%3 = 0, 


and so we can find any number of unknown quantities, if 
each equation contains one unknown quantity more than those 
before. 

In the general case when n unknown quantities are to be 
determined from n linear equations each equation will contain 
all the unknown quantities, and therefore we cannot find them 
one after the other as in the case just treated. But it can be 
shown that by means of very simple constructions the general case 
is reduced to a set of equations, such as has just been treated. 

Let us begin with two 
equations and two un- 
known quantities. 


do + 21 + ax = 0, 
Fic. 19. bo t bya, + bow, = 0. 


The lines ao, a1, a2 are laid off on a horizontal line 0494142 and 
the lines bo, 61, b: on another horizontal line O’BoB,B, (Fig. 19). 
Now let us join O and 0’, Ap and Bo, Ai and By, Ap and Bz by 
straight lines and let us draw a third horizontal line intersecting 
them in the points 0’’C,C,C,. These points correspond to a 
certain linear function 


Co + C1%1 + Coa, 


and it can be shown that it vanishes when x, and 2 are the same 
values for which the first two linear functions vanish. Let the 
distance of the first two horizontal lines be / and the distance of 


the third from the first and second h and k. Then it can readily 
be seen that 


h k h 
co = a + 7 (bo — aa) = 7 ao + 5 bo. 
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For a parallel to OO’ through Ap defines with the line AoBo on 
the third and second horizontal line segments equal to ey) — ao 


and bo — ap and as these segments have the ratio h/l, it follows 
that 


h k h 
Co = dot 7 (bo — ao) = 7 ao + 7 bo 


By drawing a parallel to 4oBy through A; and to A,B, through 
A, or through B, (which comes to the same thing), we convince 
ourselves in the same way that 


a=at ea) = Kats, 
and 
h 
= a+ 7 (b: — 2) = 
Multiplying the equation 
ao + Ut + ar = 0 
by &/l and the equation 
bo + bia + bor, = 0 
by h/l and adding the two products, we obtain 


Co + C101 + may = 0. 


The third horizontal need not lie between the first two. If it 
lies below the second we have merely to give k a negative value 
and if it lies above the first we have to give h a negative value 
and the same formule for ¢o, ¢, ¢@ hold good. Consequently the 
conclusion remains valid, that from the first two equations the 
third follows. 

Now as we are perfectly at liberty to draw the third horizontal 
line where we please, we can let it run through the intersection 
of the straight lines 4:B; and A2B:. In this case the points Ci 
and C; must coincide and consequently c must vanish. If q 
does not vanish we can by what, has been shown above find 7 
and with 2, we can find x from either of the two first horizontal 
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lines. In case c also vanishes, that is to say, in case the three 
straight lines A.B2, A1B;, AoBo all pass through the same point, 
while OO’ does not pass through it, the two given equations 
cannot simultaneously be satisfied. For if they were, it would 
follow that 

Cot 1% + Gar = 0, 


and as c, and @ are zero ¢) would have to be zero, which it is not 
as OO’ is supposed not to pass through the intersection of ABs, 
A,B; and AopBo. If on the other hand all four lines A.B, AiBi, 
AoBo, OO’ pass through the same point, ¢o, ¢; and ¢ will all three 
vanish. In this case the two given equations do not contradict 
one another, but bob:b. will be proportional to aoaia2. The 
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second equation will therefore contain the same relation between 
x, and 2 as the first, so that there is only one condition for 2, 
and 2 to be satisfied. We may then assign any arbitrary value 
to one of them and determine the value of the other to satisfy the 
equation. 

In the case of two linear equations of any number of quantities 
X1, Ve, *** Ln we can by the same graphical method eliminate one 
of the quantities. In Fig. 20 this is shown for two linear equa- 
tions with six unknown quantities. The two horizontal lines 
OApA1A2A3A4A5A6 and O' BoB, B,B;B1B;B, represent two linear 
equations. Through the intersection of 4;B; and A,B, a third 
horizontal line is drawn intersecting the lines 00’, A pBo, A1Bi, 

. ABs in OCC, «++ Cs. As Cs and C; coincide, the line c, 
vanishes and 2 is eliminated, so that the equation assumes the 
form 


Co + C121 PF de eg%3 + Cots + Cote = 0. 


i 


> 


, 


— 
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Suppose now that a set of six equations with six unknown quan- 
tities is represented geometrically on six horizontal lines. We shall 
keep one of these; but instead of the other five we construct five 
new ones from which one of the unknown quantities has been 
eliminated by means of the first equation. Now it may happen 
that at the same time another unknown quantity is eliminated, 
then this quantity remains arbitrary. Of the five new equations 
we again keep one that contains another unknown quantity and 
replace the four others again by four new ones from which this 


unknown quantity has been eliminated. Going on in this 


manner the general rule will be that with each step only one 
quantity is eliminated, so that at last one equation with one un- 
known quantity remains. Instead of the given six equations 
with six unknown quantities each, we now have one with six, 
one with five and so on down to one with one. The geometrical 
construction shows that this system is equ'valent to the given 
system, for we can just as well pass back again to the given 
system. We have seen above how the unknown quantities 
may now be found geometrically. It may however happen in 
special cases that with the elimination of one unknown quantity 
another is eliminated at the same time. To this we may then 
assign an arbitrary value without interfering with the possibility 
of the solution. Finally all unknown quantities may be elimi- 
nated from an equation. If in this case there remains a term 
different from zero it shows that it is impossible to satisfy the 
given equations simultaneously. If no term remains, the two 
equations from which the elimination takes its origin contain the 
same relation between the unknown quantities and one of them 
may be ignored. 

§ 5. The Graphical Handling of Complex Numbers.—A complex 
number 

z=a+y 


is represented graphically by a point Z whose rectangular coordi- 
nates correspond to the numbers a and y. The units by which 
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the codrdinates are measured, we assume to be of equal length. 

We might also say that a complex number is nothing but an 

algebraical form of writing down the codrdinates of a point ina 

plane. And the calculations with complex numbers stand for 

certain geometrical operations with the points which correspond 

to them. 
By the “sum” of two complex numbers 


a=uatyt and » = 2+ yt 


we understand the complex number ; 


%3 = x3 + Yysd 
where 


m= 1+a and y= y+ y%, 
and we write 


2 = 2+ &. 


Graphically we obtain the point Z3; representing z; from the 
points Z, and Z, representing 2, and x by drawing a parallel 
to OZ, through Z; and making Z,;P 
(Fig. 21) equal to OZ, in length 
and direction or by drawing a paral- 
lel through Z, and making Z,P 
equal to OZ; in length and direc- 
tion. The coordinates of P are 
evidently equal to a+ 2 and 
Yi + gy. 

Two complex numbers z and 2’ 


Fig. 21, are called opposite, when their sum 
is zero. 
z+2)=0 or e=—-2' and y= —y' or g=—2. 


The corresponding points Z and Z’ are at the same distance from 
the origin O but in opposite directions. 

The difference of two complex numbers is that complex 
number, which added to the subtrahend gives the minuend 


2 + (2 — &) = 2. 


ie i 


Lae 
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Therefore 
Zs — 2% = (m — %) + (Hm — Hm). 
This may also be written 


Z1 + 2’ where »’= —m% = — 2m — yo. . 


That is to say, the subtraction of the complex number x from x 
may be effected by adding the opposite number — m. For the 
geometrical construction of the point Z corresponding to 2, — x 
we have to draw a parallel to OZ, through Z, and from Z, in 
the direction from Z, to O we have to lay off the distance Z.0. 
Or we may also draw from 0 a line equal in direction and in length 
to Z,Z;. ‘This will also lead to the point Z representing the 
difference 21 — %. 

The rules for multiplication and division of complex numbers 
are best stated by introducing polar codrdinates. Let r be the 
positive number measuring the distance OZ in the same unit 
of length in which 2 and y measure the abscissa and ordinate, so 
that 

r=V2+% 


and let ¢ be the angle between OZ and the axis of x, counted in 
the direction from the positive axis of x toward the positive 
axis of y through the entire 
circumference (Fig. 22). Then 
we have 

z=rceosg, y=rsing 


and 


z=a+yi=r(cose+ sin gi). 


Let us call r the modulus 
and ¢ the angleof z. The an- Fic. 22. 
gle may be increased or di- 
minished by any multiple of four right angles without altering 
z, but any alteration of r necessarily implies an alteration of z. 
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According to Moivre’s theorem, we can write 
a= re™ 
By the product of two complex numbers 
a= ne’ and x» = ners 


we understand that complex number 2; whose modulus 13 is 
equal to the product of the moduli r; and r, and whose angle ¢s 
is the sum of the angles ¢; and ~ or differs from the sum only by 
a multiple of four right angles 


23 = Am = Tyree Pit 24, 


The definition of division follows from that of multiplica- 
tion. The quotient 2 divided by ~% is that complex number, 
which multiplied by ~ gives x. Therefore the product of its 
modulus with the modulus of z must be equal to the modulus of 
2, and the sum of its angle with the angle of x must be equal to 
the angle of x. Or we may also say the modulus of the quotient 
2/2 is equal to the quotient of the moduli r;/rz and its angle is 
equal to the difference of the angles g, — ~. An addition or 
subtraction of a multiple of four right angles we shall leave out 
of consideration as it does not affect the complex number nor 
the point representing it. 

The geometrical construction corresponding to the multi- 
plication and division of complex numbers is best described by 


considering two quotients each of two complex numbers that 
give the same result. Let us write 


21/2 = 23/24. 


The geometrical meaning of this is that 


1/T2 = 13/T4, 
and 


Pl ime Cs a OS eee Od. 


That is to say, the triangles Z,0Z, and Z;0Z, are geometrically 


3 


L 
? 
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similar (Fig. 23). When three of the points Z;, Z2, Z3, Zs are 
given the fourth can evidently be found. For instance let 
Z;, Z2, Zs be given. Draw a parallel to Z,Z, intersecting OZ, 
at a distance 7, from O. This point together with the inter- 
section on OZ; and with O will form the three corners of a tri- 
angle congruent to the triangle Z,Z;0. It will be brought into 


Z4=1 


ic. 23. Fie. 24. 


the position of Z,;Z30 by being turned round 0 so as to bring the 
direction of the side in OZ, into the position of OZ, Thus the 
direction of OZ; and its length may be found. 

This construction contains multiplication as well as division as 
special cases. Let Z, coincide with the point x = 1, y = 0, so 
that z, = 1 (Fig. 24), then we have 


21/20 = 23 OF YB = zs. 


From any two of the points Z, Z2, Z; a simple construction gives 
us the third. 
The geometrical representation of complex numbers may be used 
to advantage to show the properties of harmonic oscillations. 
Let a point P move on the axis of 2, so that its abscissa at the 
time tis given by the formula 


x=r cos (nt + a), 


n,rand a being constants. We call r the amplitude and nf + a 
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the phase of the motion. The point P moves backwards and 
forwards between the limits =r and a=—r. The time 
T = 2n/n is called the period of the oscillation, it is the time in 
which one complete oscillation backwards and forwards is per- 
formed. 
Now instead of a let us consider the complex number 
z= reos(nt+ a) + rsin(nt+ a) 


or ; 
z= rete 


of which z is the abscissa and let us follow the movement of 
the point Z. For ¢ = 0 we have 


z= rem, 
Designating this value by zo, we can write 
3 =" 2e™. 
The geometrical meaning of the product 
zoenti 


is that the line OZ) is turned round O through the angle nt. For 
the modulus of e”” being equal to 1 the modulus of 29 is not 
changed by the multiplication. The 
movement of the point Z therefore 
consists in a uniform revolution of 
OZ round 0. At the moment t= 0 
the position is OZ) and after the 
time 7’ = 2z/n the same position is 
occupied again. The revolution goes 
on in the direction from the positive 
axis of a to the positive axis of y 

(Fig. 25). 
Sr FS The movement of Z is evidently 
simpler than the movement of the 

projection P of Z on the axis of 2. 

Let us consider a motion composed of the sum of two harmonic 
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motions of the same period but of different amplitudes and 
phases 
x = 7, cos (nt + ay) + 7 cos (nt + a), 


and let us again substitute the motion of the point Z correspond- 
ing to the complex number 


z= rete 1 ppglnttands, 


For ¢t = 0 the first term is 


2 = res 


and the second term 
m= = ne, 
Introducing 2 and x into the expression for z we have 
z= gent — gent? — (2 + Bent —_ agent 
where 
™ r 243 = 21 + 226 
This shows at once that the movement of Z is a uniform circular 
movement consisting in a uniform revolution of OZ round 0. 
The position at the moment t = 0 is OZ; corresponding to the 


complex number 
23 = AB + 29. 


The projection of Z on the axis of x has the abscissa 
x = 73 cos (nt + a3) 


where 7; and a; designate modulus and angle of z. Thus the 
sum of two harmonic motions of the same period is shown also 
to form a harmonic motion. 

The same holds for a sum of any number of harmonic motions 
of the same period. For the complex number 


g = nemtteni 4 ppelnttadi t ... 4 pelnttays 
where 71, 72, °°* Tj} Q1, @,-°+* a and m are constants may be 
written : 

a zen -L zen a. ‘apa a (sed 


or 
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f= set, 
where 


Z=Aateate--+% 


The movement of Z therefore, excepting the case z9 = 0, consists 
in a uniform revolution of OZ round O, OZ always keeping the 
same length equal to the modulus of zo. The position of OZ at 
the moment ¢ = 0 is OZp. 

The motion of a point P whose abscissa is 


x = ae cos (nt + a) 


where a, k, n, a are constants (a and k positive) is called a damped 
harmonic motion. It may be looked upon as a harmonic motion, 
whose amplitude is decreasing. To study this motion let us 
again substitute a complex number 


z = ae cos (nt + a) + ae™ sin (nt + a)i, 
or 
2 = aeH.elattadé, 


or 
— Boe se 


where 2p is written for the complex constant ae™. 


The product 
ze 


is a complex number corresponding to a point Z; on the same 
radius as Zo, coincident with Zp at the moment ¢ = 0 but ap- 
proaching O in a geometrical ratio after¢ = 0. In unit of time 
the distance of Z; from O decreases in the constant ratio e~* : 1. 
The multiplication with e™* turns OZ; round O through an angle 
nt. We may therefore describe the motion of Z as a uniform 
revolution of OZ round O, Z at the same time approaching O 
at a rate uniform in this sense that in equal times the distance 
is reduced in equal proportions (Fig. 26). At the moment ¢t = 0 
the position coincides with Zp. We speak of a period of this 
motion meaning the time 7 = 2z2/n in which OZ performs an 
entire revolution round OQ, although it does not come back to its 


~ 
oe 
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original position. Any part of the spiral curve described by Z _ 


in a given time is geometrically similar to any other part of the 


curve described in an interval of equal duration. For suppose 
the second interval of time hap- 
pens 7 units of time later, we 
shall have for the first interval 


z= ze Bt. enté, 
and for the second interval 
2’ = zoe * tt) «en (ttt) 
Now if 2: and z are the values 
of z at two moments 4, and & of 
the first interval and x’ and ’ Fig. an 


the corresponding values of 2’ 
at the moments 4 + 7 and & + 7 of the second interval, we have 


, 
21 — g—k(ti-—ts) . en (a—te)t — all 


Ae 
a) 2) 


Therefore the triangle Z:;0Z, is geometrically similar to the 
triangle Z;/OZ,’. As Z, and Z, may coincide with any points 
of the first part of the curve, the two parts are evidently geo- 
metrically similar. 

The projection of Z on the axis of 2 performs oscillations 
decreasing in amplitude. The turning-points correspond to those 
points of the spiral curve described by Z, where its tangent is 
parallel to the axis of y, that is to say, where the abscissa of dz/dt 
vanishes. 


Now 
e = 2o(— k+ nije Me™ = (—k + mi)z 
or 
dz 
dF yt 


34 GRAPHICAL METHODS. 


where p and ) are the modulus and angle of the complex number 
—k+ ni. 

Consequently, if we represent dz/dt by a point Z’, the triangle - 
Z'OZ will remain geometrically similar to itself. The turning 
points of the damped oscillations correspond to the moments 
when OZ’ is directed vertically upward or downward or when the 
angle of dz/dt is equal to 7/2 or 37/2. The angle of z will then 
be 7/2 — d or 37/2 — X plus or minus any multiple of 27. As 
the angle of z, on the other hand, is changing in time according 
to the formula 

nt + a, 


we find the moments where the movement turns by the equation 


nt+a= 7r/2—r+ 2Nz, 
or 
nt + a = 32/2 —A+ 2Nz, 


N denoting any positive or negative integral number. The time 
between two consecutive turnings is therefore equal to z/n, that 
is, equal to half a period. All the points Z corresponding to 
turning points lie on the same straight line through the origin O 
forming an angle 37/2 — \ with the direction of the positive axis 
of a. The amplitudes of the consecutive oscillations therefore 
decrease in the same proportion as the modulus of 2, that is 


to say, in half a period in the ratio ee 

Let us consider the vibrations of a system possessing one 
degree of freedom when the system is subjected to a force varying 
as a harmonic function of the time and let us limit our considera- 
tions to positions in the immediate neighborhood of a position 
of stable equilibrium. If the quantity z determines the position 
of the system the oscillations satisfy a differential equation of 
the form 


x dx 
m 5 + ka + nx = F cos (pt)! 


where m, k, n, p, F are positive constants. 


1 See for instance Rayleigh, Theory of Sound, Vol. I, chap. III, § 46. 


E 
. 
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This is another case where the introduction of a complex 


z=aty 


and the geometrical representation of complex numbers helps to 
form the solution and to survey the variety of phenomena that 
may be produced. 

‘In order to introduce z let us simultaneously consider the 
differential equation 


a d 
moat ket ny = F sin pt, 

and let us multiply the second equation by 7 and add it to the 
first. We then have 


Pe ee pines 
mat kat wa = Frere. 


The movement of the point Z representing the complex number 
z then serves as well to show the movement corresponding to x. 
We need only consider the projection of Z on the axis of a. 

A solution of the differential equation may be obtained by 
writing 

z= 290", 
Introducing this expression for z and cancelling the factor e?” 
we have 
zo(— mp? + kpi + n*) = F, 
F 
== mp? + kpi+ n? 


or 


zo is a complex constant, that may be represented geometrically 
as we shall see later on. 


This solution 
ye) ze?” 


is not general. If z’ denotes any other solution so that 
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mo hE nite = Pent 


we find by subtracting the two equations 


noe?) db pe “ — 2) 


or writing 


+ n?(z’ —2z) =0 


The general solution of this equation is 
“= me uv — me, 


where uw and wz are arbitrary constants and ; and 2d; are the 
roots of the equation for \ 


m2 + ky + n? = 0, 
ne a 


— + r- 
2m 4m 


If k?/4m? is greater than n, so that the square root has a real 
value, V k?/4m? — n? will certainly be smaller than k/2m. There- 
fore \; and 2 will both be negative and the moduli of the complex 
numbers we and we’ will in time become insignificant. If, 
on the other hand, k?/4m? is smaller than n?, both complex 
numbers we” and we correspond to points describing spirals that 
approach the origin, as we have seen above, in a constant ratio 
for equal intervals of time. Therefore they will also in time 
become insignificant. 
After a certain lapse of time the expression 


2 = 2 er 


will therefore suffice to represent the solution. 

The point Z moves uniformly in a circle round O of a radius 
equal to the modulus of zo, completing one revolution in the 
period 2z/p, the period of the force acting on the system. The 


eT 
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movement of the projection of Z on the axis of x is given by 
x = rocos (pt+ a), 


where ro is the modulus and a@ the angle of zo. It is a harmonic 
movement with the same period as that of the force: F cos pt, 
but with a certain difference of phase and a certain amplitude 
depending on the values of F, m, k, n, p. 

It is important to study this relation in order to survey the 
phenomena that may be produced. For this purpose the geo- 
metrical representation of complex numbers readily lends Hey 

In the expression for 2 

F 
— mp +kpi+ n” 


20 = 


let us consider the denominator 
— mp* + kp + n’, 


and let us suppose the period of the force acting on the system 
not determined, while the constants of the system m, k, n and 
the amplitude of the force F have given values. The quantity p 
is the number of oscillations of the force during an interval of 
2r units of time. This quantity p we suppose to be indeter- 
minate and we intend to show how the amplitude and phase 
of the forced vibrations compare with the amplitude and phase 
of the force for different values of p. 

Let us plot the curve of the points corresponding to the complex 


number 
n? — mp? + kpr, 


where p assumes the values p = 0 to + ©. 

This curve is a parabola whose axis coincides with the axis of 
x and whose vertex is in the point z = n?, y= 0. We find its 
equation by eliminating p from the equations 


z=n—mp, y= kp, 
V1zZ., 
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™ 
ers — 4/2 
0 — FL ye" 


But it is better not to eliminate p and to plot the different points 
for different values of p. In Fig. 27 the curve is drawn for p = 0 
to 3 and the points for 
p=0,1, 2, 3 are marked. 
The ordinates increase in 
proportion to p; they are 
equal to 0, k, 2k, 3k for 
o=0, 1, 2, 3. The die 
tance between the projec- 
tion of any point of the 
curve on the axis of x and 
the vertex is proportional to p*. It is equal to 0, m, 4m, 9m for 
n= 0,122, 3, 

For any point P on the parabola let us denote the distance 
from O by r and the angle between OP and the positive axis of 
x by ¢ so that 


n? — mp? + kpi = re®*. 
Then we have 
Zp. Sore", 
and consequently 


F 
z= peas 
and 


F 
«=~ cos (pt — ¢). 


The amplitude F/r of the forced vibration is inversely propor- 
tional to r. Thus our Fig. 27 shows us what the period of the 
force must be to make the forced vibrations as large as possible. 
It corresponds to the point on the parabola whose distance from 
O is smallest. It is the point where a circle round O touches the 
parabola. In Fig. 27 this point is marked R. It may be called 
the point of maximum resonance. When the constants of the 
system are such that the ordinate of the point, where the parabola 
intersects the axis of y is small in comparison with the abscissa 


a 


. quarter of a period. Keep- 
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of the vertex, then OR will lie close to the axis of y (Fig. 28). In 
this case the angle between OR and the positive axis of x will be 
very nearly equal to 90°, that is to say, the forced oscillations will 
lag behind the force oscil- 
lations by a little less than a 


ing m and n constant, this 
will take place for small val- 
ues of k, 7@. e., for a small 
damping influence. A small Fie. 28. 
deviation of p from the fre- 
quency of maximum resonance will throw the point P awayfrom R, 
so that r increases considerably and ¢ becomes either very small 
(for values of p smaller than the frequency of maximum resonance) 
or nearly equal to 180° (for values of p larger than the frequency 
of maximum resonance). In other words for small values of k the 
maximum of resonance is very sharp. A deviation of the period 
of the force from the period of maximum resonance will lessen the 
amplitude of the forced vibration considerably. The lag of its 
phase behind that of the force will at the same time nearly vanish, 
when the frequency of the force is decreased or it will become nearly 
as large as half a period, when the frequency of the force is in- 
creased. For larger values of k the parabola opens out and this 
phenomenon becomes less marked. The minimum of the radius r 
becomes less pronounced. ‘The angle between OR and the axis of 
x becomes smaller and smaller and for a certain value of & and all 
larger values the point R will coincide with the vertex of the para- 
bola.. In this case, there is no resonance. When the period of 
the force increases indefinitely (p becoming smaller and smaller) 
the amplitude of the forced vibration will increase and will 
approach more and more to the limit 

F 


nm? 


but there will be no definite period for which the forced vibra- 
tions are stronger than for all others. 


CHAPTER II. 


Tue GRAPHICAL REPRESENTATION OF FUNCTIONS OF ONE OR 
Morr INDEPENDENT VARIABLES. 


§ 6. Functions of One Independent Variable.—A function y of 


one variable x 


y = fe) 


is usually represented geometrically by a curve, in such a way 
that the rectangular codrdinates of its points measured in certain 
chosen units of length are equal to 2 and y. This graphical rep- 
resentation of a function is exceedingly valuable. But there is 
another way not less valuable for certain purposes, more used in 
applied than in theoretical mathematics, which here will occupy 
our attention. 

Suppose the values of y are calculated for certain equidistant 
values of 2, for instance: 
x= — 6,—5, —4, — 3, —2, — 1,0, 

+ 1, -+ 2, 4:3, +> 4) ee 
and let us plot these values of y in a uniform scale on a straight 
line. Draw the uniform scale on 
one side of the straight line and 
mark the points that correspond 
to the calculated values of y on 
the other side of the straight line. 
Denote them by the numbers z 
that belong to them (Fig. 29). 
The drawing will then allow us to 
read off the value of y for any of 

Pie29. the values of 2 with a certain ac- 

curacy depending on the size of the 

scale and the number of its partitions and naturally on the fine- 
40 
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_ may be imagined to signify the 
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ness of the drawing. It will also allow us to read off the value 
of y for a value of x between those that have been marked, if 
the intervals between two consecutive values of a are so small 
that the corresponding intervals of y are nearly equal. We can 
with a certain accuracy interpolate values of 2 by sight. On the 
other hand, we can also read off the values of a for any of the 
values of y. We shall call this the representation of a function 
by a scale. 

We can easily pass over to the representation of the same 
function by a curve. We need only draw lines perpendicular 
to the line carrying the scales through the points marked with 
the values of x and make their length measured in any given 
unit equal to the numbers 2 that correspond to them (Fig. 29). 

In a similar way we can pass 
from the representation of the 
function .by a curve to the rep- 
resentation by a scale. 

The representation by a scale 


Position 


movement of a point on a straight 
line, the values of x meaning the 
time and the points marked with 
these values being the positions 
of the moving point at the times 
marked. By passing over to the curve the movement in the 
straight line is drawn out into a curve with the time as abscissa 
(Fig. 30). 

The representation by a scale is used in connection with the 
representation by a curve for the purpose of drawing a function 
of a function. 

Let y be a function of x and 2 a function of t. Then we wish 
to represent y as a function of ¢. 

Let y = f(x) be given by a curve in the usual way and let 
x = ¢(t) be given by a scale on the axis of 2 marking the points 
where t = 0, 1, 2, ---, 12. We then find the values of y corre- 
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sponding to the values t = 0, 1, 2, --+, 12 by drawing the ordi- 
nates of the curve y = f(a) for the abscissas marked t = 0, 1, 2, 


—---y =f (¢(t)) 
—-y=f (@) 


-++, 12. These ordinates as a rule will not be equidistant. But 
as soon as we move them so as to make them equidistant, they 
form the ordinates of the curve 


y = f(e@®) 


with ¢ as abscissa (Fig. 31). 

The representation of a func- 
tion by a scale may be general- 
ized in the respect that neither 
of the two scales facing one an- 
other on the straight line need 
necessarily be uniform. The in- 
tervals of both scales may vary 
from one side of the scale to the 
other. If the variation is suffi- 
ciently slow the interpolation can nevertheless be effected with 


accuracy. We may look at this case as composed of two cases 
of the first kind. 


f(z) =y and y=. 


a 


lt Ri 
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These scales are placed together, so that the scale x touches the 


scale ¢ 


F(x) = 9, 


while the scale y is cut out (Fig. 32). 

$7. The Principle of the Slide Rule.—Let us investigate how 
the relation between x and ¢ changes by sliding the a- and ¢-scales 
along one another. 

If we slide the x-scale through an amount y = ¢ so that a 
point of the x-scale that was opposite to a certain point y of the 
y-scale, now is opposite y + c, then the relation between zx and t 
represented by the new position of the scales will be given by 
the equation 


f@) =9@ +e. 


If x, ¢ and 2’, t’, denote two pairs of values that are placed 
opposite to one another, we shall have simultaneously 


f(z) =g +e, 
fx’) = 9@/) +e, 
or by eliminating c 


~ f@) — 9 =f) — gf). 


The ordinary slide rule carries two identical scales y = log x and 
y = log ¢ that are able to slide along one another, 2 and ¢ running 
through the values 1 to 100. We therefore have 


log x — logt = log a’ — log?’, 


or 


That is to say, in any position of the 2- and ¢-scale any two values 
zx and t opposite each other have the same ratio (Fig. 33), This 
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is the principle on which the use of the slide rule is founded. 


It enables us to calculate any of the four quantities a, t, 2’, t’. 


if the other three are given. Suppose, for example, 2, ¢, 2’ 
known. We set the scales so that x appears opposite to 1, 


then ?’ is read off opposite to 2’. On the other edges the slide 


rule carries two similar scales one double the size of the other 


(Fig. 34). We may write 
y=2logX and y= 2 log T. 


By means of a little frame carrying a crossline and sliding over 
the instrument, we can bring the scales x and T or t and X op- 
posite each other. If, for example, for any position of the 
instrument 2, 7 and x’, T’ are two pairs of values opposite each 
other, then 

log x — 2 log T = log 2’ — 2 log T’, 
or 


=> 7? 4 


ale 


If any three of the four quantities z, 7, 2’, T’ are known the 
fourth may be read off. Thus we find the value 


eT” 
Te 


by setting 7’ opposite to x and reading off the value opposite to 
T’. Or we can find the value of 


Ja’ 
Var 
; a 
by setting x opposite to T and reading off the value opposite 2’. 


: 
] 


\ 


. 


| 


/ 
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Let us reverse the part that carries the scales t, 7 so that x 
slides along 7 and X along ¢, but in the opposite order (Fig. 35). 


Fie. 35. - 
The scales ¢, T may then be expressed by 
y=Il—logit and y=1— 2logT, 


l being the entire length of the scales. 

By setting the instrument to any position and considering the 
scales x and ¢ or X and T by means of the cross line we have 
logx + logt = logz’ + log?’ and log X + log T = log X’ + log T’ 
or 

a= and XT = X'T’, 


so that any two values opposite to one another have the same 
product. 
For x and T we have 


log 2 + 2log T = log a’ + 2 log 7’, 
or 
aT? = 2/T” 


Let us apply this to find the root of an equation of the form 


w+ au= b. 

Divide by wu so that : 
2 eape 

w+a if 


and set T = 1 opposite to X = 6. Then taking T = uw we find 
on the same cross line tf = wv? and X = b/u, so that we read the 
two values wu? and b/w directly opposite to each other on the 
scales t and X. If b/w is positive, it decreases while wv? increases. 
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Running our eye along we have to find the place where the differ- 
ence b/u — u? is equal to a. Having found it the T-scale gives 
us the root of the equation. For example take 


uw — 5u = 3, 
or 


We set T = 1 opposite XY = 3 and run our eye along the scales 
X and t (Fig. 36), to find the place where t — 5 = X. We find 


6 cT 
orese29S F & F 


OF OF 06 


it approximately at t = 6.2, and on the T-scale we read off 
T = 2.50 as the approximate value of the root. This is the 
only positive root. But for a negative root 3/u is negative, 
and therefore the positive value of 3/w plus u? would have to be 
equal to 5. We run our eye along and find t = 3.37 opposite to 
X = 1.63, approximately corresponding to JT = 1.84. There- 
fore — 1.84 is another root. As the coefficient of wu? in the first 
form of the equations vanishes it follows that the sum of the 
three roots must be equal to zero. This demands a second 
negative root approximately equal to — 0.66. To make sure 
that it is so, we set the instrument back and take the other end 
of the T-scale as representing the value J = 1 and give it the 
position this end had before. Running our eye along the 
scales X and t, we find ¢ = 0.43 opposite to X = 4.57, giving 
X +t= 5.00. On the T-scale we find 0.655, so that the third 
root is found equal to — 0.655. 

When 6 is negative there is always one and only one negative 
root. For w running through the values wu = 0 to — 00, w2—b/u 
will run from — 00 to + oo without turning. When 6 is positive 
there is always one and only one positive root; for then w2 — b/u 
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runs from — © to + oo for w= 0 to +o, In the first case 
there may be two positive roots or none; in the second case there 
may be two negative roots or none. For positive values of a 
one root only exists in either case. This is easily seen in the first 


form of the equation 
wu + au = b, 


because from a positive value of a it follows that wu? + aw will 
for wu = — oo to + ©, run from — o to + oo without turning 
and will therefore pass any given value once only. 

In order to decide whether in the case of a negative value of a 
there are three roots or only one let us write 


For negative values of b we have to investigate whether there 
are positive roots. For positive values of wu the function w?—b/u 
has a minimum, when the differential coefficient vanishes, 7. e., for 


ee , 
uw 


or 


Having set our slide rule so that ¢ gives us uw? and X gives us 
— b/u, we find the value w where the minimum takes place by 
running our eye along and looking for the values X, ¢ opposite 
each other for which X is twice the value of t 


2t = X. 


Then ¢-++ X is the minimum of wu? — b/u, so that there will be 
two or no positive roots according to t+ X being smaller or 
larger than — a. For positive values of 6, we have to find out 
whether there are negative roots. The criterion is the same. 
After having set J’ = 1 opposite to b and having found the 
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positive root, we find the place where 
2t = X. 


Then ¢ + X is the minimum of all values that u? — 4/u assumes 
for negative values of wu. If the minimum is smaller than — a 
there are two negative roots; if it is larger there are none. If it 
is equal to — a the two negative roots coincide. 

For the equation 


Lilw 


for instance, we find ¢t = 1.31 opposite to X = 2.62 (Fig. 36), 
so that 2¢ = 2.62 = X. Now t+ X = 3.93 is smaller than 5, 
therefore wu? — 3/u will assume the value 5 for two negative 
values of uw on either side of the value w= — T = — 1.143 
for which the minimum of wu? — 3/w takes place. 

On the same principle as the slide rule many other instruments 
may be constructed for various calculations. In all these cases 
we have for any position of the instrument 


f(z) — 9® = f(z’) — 9), 


where x, ¢ are any readings of the two scales opposite each other 
and x’t’ the readings at any other place. f(x) and g(t) may be 
any functions of x and ¢. It will only be desirable that they 
be limited to intervals of x and ¢, which contain no turning 
points. Else the same point of the scale corresponds to more 
than one value of x or ¢ and that will prevent a rapid reading 
of the instrument. 

Let us design an instrument for the calculation of the increase 
of capital at compound interest at a percentage from 2 per cent. 
upward. If 2 is the number of per cent. and ¢ the number of 
years, the increase of capital at compound interest is in the pro- 


portion 
x t 
(1+ 35). 


Ee ee 
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We can evidently build an instrument for which 


(8 io) = (1+ ia) 


For taking first the logarithm and then the logarithm of the 
logarithm, we obtain 


log ¢ + log log (1 = Te log t’ + log log (1 ao i sa) 


We have only to make the z-scale 


y = + log log (1 + a — log log (1 + ma) 


and the ¢-scale 
y = log n — logt. 


For « = 2 we have y = 0 and therefore in the normal position 
of the instrument ¢= nm. On the other end we have ¢ = 1 and 
therefore y = logn. Now let us take n = 100, so that y = 2 
for t= 1. Say the length of the instrument is to be about 24 
em., then the unit of length for the y-scale would have to be 12 
em. In the normal position of the instrument the readings 2, t 
opposite to each other satisfy the equation 


t 9 100 
(14 TH (1+ aa 


Opposite t = 1, we read the value 2, = 624 and this gives us 


100 
os Se: D) 
( i ma) =e a 1+ 6.24 = 7.24. 


A capital will increase in 100 years at two per cent. compound 
interest in the proportion 7.24:1. Or we may also say the 
number 2x; = 624 read off opposite t = 1 is the amount which is 
added to a capital equal to 100 by double interest of 2 per cent. 
in 100 years. The same position of the instrument gives us the 
number of years that are wanted for the same increase of capital 


5 
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at a higher percentage. For all the values 2, ¢ opposite to each 
other satisfy the equation 


a 8 
(: + a) = 7.24. 


For any other given percentage x and any other given number 
of years ¢ the increase of capital is found by setting x opposite 
to t and reading the z-scale opposite tof = 1. The only restric- 
tion is that the ratio is not greater than 7.24, else t= 1 will 
lie beyond the end of the z-scale. 

For a given increase of capital the instrument will enable us 
either to find the number of years if the percentage is given, or 
the percentage if the number of years is given, subject only to 
the restriction mentioned. 

We can build our instrument so as to include greater increases 
of capital by choosing a larger value of n. n = 1000, for in- 
stance, will make y = 3 fort = 1. If the instrument is not to 
be increased in size the scales would have to be reduced in the 
proportion 2 : 3. 

Let us consider another instance 

1 1 


y= yn 


1 

ae 

In the normal position of the instrument the scale division 
marked x = © corresponds to y = 0 and is opposite to t = n. 
If we have ¢ = 00 on the other end, the length of the instrument 
will correspond to y = 1/n. Let us choose n = 0.1, so that the 
length of the instrument is y = 10. That is to say, the unit of 
length of the y-scale is one tenth of the length of the instrument. 
For any position of the instrument we have 


If the scale division marked 2 = © is opposite to ¢ = ¢ we can 
write 2’ = oo, t’ = c¢ and have 


4 | 
mi 

. 
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The instrument will therefore enable us to read off any one of 
the three quantities 2, t, c, if the other two are given, the only 
restriction being that all three lie within the limits 0.1 to oo. 
The instrument may be used to determine the combined resistance 
of two parallel electrical re- 
sistances, for the resistances 
satisfy the equation 


Fie. 37. 


Similarly it may be used 

to calculate the distances of an object and its image from the 
principal planes of any given system of lenses. For if f is the 
focal length and z and ¢ the distances of the object and its im- 
age from the corresponding principal planes (Fig. 37), the equa- 
tion is 

Prank 

es ba 
On the back side of the movable part of an ordinary slide rule 
there generally is a scale 


y = 2+ log sint. 


When this part is turned round and the scale is brought into 
contact with the scale 
y = loga, 


we obtain for any position of the instrument 


log x — log sint = log 2’ — log sin?’, 


or 
x x’ 


sinf sin i’’ 


for any two pairs of values x, ¢ that are opposite each other. 


bz GRAPHICAL METHODS. 


Given two sides of a triangle and the angle opposite the larger 


of the two the instrument gives at once the angle opposite the + 


other side. Similarly when two angles and one side are given, 
it gives the length of the other side. 
If x’ = ais the value opposite to t/ = 90°, we have 


z= asini, 


Thus we can read the position of any harmonic motion for any 
value of the phase. 
An instrument carrying the scales 


y=logsine and y= logsint 
enables us to find any one of four angies 2, ft, x’, t’ for which 


sinz sin 2’ 
sné sin v 


if the other three are given. Thus, knowing the declination, 
hour angle and height of a celestial body, we can read the azimuth 
on the instrument. We have only to take x = 90° — height, 
t = hour angle, 2’ = 90° — declination, then ¢/ = azimuth or 
180° — azimuth. 

It is not necessary to carry out the subtraction 90° — height and 
90° — declination. The difference may be counted on the scale 
by imagining 0° written in the place of 90°, 10° in the place of 
80° and so on and counting the partitions of the scale backwards 
instead of forward. 

§ 8. Rectangular Codrdinates with Intervals of Varying Size.— 
The two methods of representing the relation between two 

variables either by a curve connecting the codrdinates or by 
scales facing each other lead to a combination of both. 


Suppose the rectangular codrdinates x and y are functions of 


u and 2, 
x=g(u) and y= yr). 


The function x = ¢(u) is represented by a uniform scale for x 
on the axis of abscisse facing a non-uniform scale for u. The 
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function y = y¥(v) is represented by a uniform scale for y on 
the axis of ordinates facing a non-uniform scale for v. Through 
the scale-divisions wu let us draw vertical lines, and through the 
scale-divisions v let us draw horizontal lines. These two systems 
of parallel lines form a network of rectangular meshes of various 
sizes (Fig. 38), and any equation between uw and v may be TEP 
sented by a curve in this plane. 

The usefulness of this method will be seen by some examples. 
It enables us by a clever choice of the functions g(u) and y(v) 


Fig. 38. 


to simplify the form of the curve. It is easily seen, for instance, 
that a curve representing an equation f(u, v) = 0 may always be 
replaced by a straight line, if we choose the w-scale properly. 
For when the points u = 1, 2, 3, 4, --- of the curve are not on 
a straight line, let them be moved to a straight line without 
altering their ordinates (Fig. 39). This will change the u-scale 
but it will not alter the equation f(u, v) = 0 now represented by 
the straight line. 
Suppose we want to represent the relation 


where a and b are given numbers. If w and v were ordinary 
rectangular coérdinates the curve would be an ellipse. But if we 


make 
z=w? and y=v 
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the equation of the line in rectangular coérdinates becomes 
v y 
atyo) 


and the curve will therefore be a straight line running from a — 
point on the positive axis of 2 to a point on the positive axis of 
y. The point on the axis of x corresponds to-the value uw = + a 
on the w-scale, and the 
point on the axis of y cor- 
responds to the value » = 
+ b on thev-scale (Fig. 40). 

Any point on the straight 
line corresponds to four 
combinations + wu, +v;—u, 
+ v; u,—v; —Uu, — 2, be- 
cause x has the same values 
for opposite values of wu 
and y for opposite values of v. We can read v as a function of 
u or u as a function of v. 

If a second equation 


Fie. 40. 


is given, we find the common solutions of the two equations by 
the intersection of the corresponding straight lines. Fig. 40 
shows the solutions of the two equations 


Uw vy 
pty a! 
and 
u 2242 
ere 


approximately equal tow = + 1.2 andv= = 2.4, 
Another function much used in mathematical physics 


cua 
Vv = Ae m2 


ee pe 
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may also be represented by a straight line by means of the same 
device. 
By making 
y=log»y, «= wv, 
we obtain 


et 
y=loga——, 


where logv and loga are the natural logarithms of v and a. 
The w-scale is laid off on the axis of x and the v-scale on the axis 
of y.and we have to join the points u = 0, v= a and u = m, 
»=a/e. The point v = a/e is found by laying off the distance 
»=1 to v=e from v= a downward (Fig. 41). We are not 
obliged to take the same units of length for z and y. 


a (U) 


Fig. 41. 


Suppose we had to find the constants a and m from two equa- 
tions 


and 


Our diagram would furnish two points corresponding to wm, %1 
and ti, %. The straight line joining these two points intersects 
the axis of ordinates at » = a and intersects the parallel through 
» = a/e to the axis of abscissee at wu = m. 
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In applied mathematics the problem would as a rule present 
itself in such a form that more than two pairs of values u, v 
would be given but all of them affected with errors of observation. 
The way to proceed would then be to plot the corresponding 
points and to draw a straight line through the points as best we 
can. A black thread stretched over the drawing may be used to 
advantage to find a straight line passing as close to the points 
as possible (Fig. 42). 

In several other cases the variables u and v are connected with 
the rectangular codrdinates x and y by the functions 


x=logu and y= log». 
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“Logarithmic paper” prepared with parallel lines for equidistant 
values of w and lines perpendicular to these for equidistant values 
of v is manufactured commercially (Fig. 43). 

By this device diagrams representing the relation 


u'r = ¢, 
where r, s, c are constants are given by straight lines. For by 
taking the logarithm we obtain | 
re + sy = loge. 
The straight ae connects the point wu = c"” on the u-scale with 
the point » = ce” on the v-scale. 
Logarithmic paper is further used to advantage in all those 
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cases where a variety of relations between the variables w and v 
are considered that differ only in w and v being changed in some 
constant proportion. If w and » were plotted as rectangular 
coérdinates the curves representing the different relations be- 
tween w and v might all be generated from one of them by altering 
the scale of the abscisse and independently the scale of the ordi- 
nates, so that the appearance of all these curves would be very 
different. Let us write 


Flu, %) = 0, 


as the equation of one of the curves. The equations of all the 
rest may then be written 


(2 5)=0 


where a, 6 are any positive constants. The points wu, v of the 
first curve lead to the points on one of the other curves by taking 
ua times as great and vb timesas great. For if we write u’ = au 
and v’ = bu the equation f(u, v) = 0 leads to the equation 


between w’ and v’: 
uo! 
f (=, q = 0. 


Using logarithmic paper the diagram of all these curves be- 
comes very much simpler. The equation f(u, v) = 0is equivalent 
to a certain equation g(x, y) = 0, where x = log u, y = log »v. 
Now let 2’, y’ be the rectangular coédrdinates corresponding to 
u’, v’ so that 


x’ = log u’ = log u+ loga = x + loga, 

y’ = log vw’ = log v+ logb = y+ log b. 
The point 2’, y’ is reached from the point 2, y by advancing 
through a fixed distance log a in the direction of the axis of x 


and a fixed distance log b in the direction of the axis of y. The 
whole curve 


flu, v7) = 90 
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drawn on logarithmic paper is therefore identical with all the 
curves 


It can be made to coincide with any one of the curves by 
moving it along the directions of z and y. 

§ 9. Functions of Two Independent Variables.—When a func- 
tion of one variable y = f(z) is represented by a curve, the values 
of x are laid off on the axis of z and the values of y are represented 
by lines perpendicular to the axis of z In a similar way a 
function of two independent variables 


a= f(z, y) 

may be represented by plotting z and y as rectangular codrdinates 
and erecting lines perpendicular to the zy plane, in all the 
points x, y, where f(z, y) is defined and making the lengths of 
the perpendiculars proportional to z. In this way the function 
corresponds to a surface in space. Now there are practical 
difficulties in working with surfaces in space and therefore it 
appears desirable to use other methods, that enable us to represent 
functions of two independent variables on a plane. This may 
be done in the following way. 

Taking 2, y as rectangular codrdinates all the points for which 
f(x, y) has the same value form a curve in the zy plane. Let 
us suppose a number of these curves drawn and marked with the 
value of f(x, y). If the different values of f(x, y) are chosen 
sufficiently close, so that the curves lie sufficiently close in the 
part of the zy plane that our drawing comprises, we are not only 
able to state the value of f(x, y) at any point on one of the drawn 
curves, but we are also able to interpolate with a certain degree 
of accuracy the value of f(z, y) at a point between two of the 
curves. As a rule it will be convenient to choose equidistant 
values of f(x, y) to facilitate the interpolation of the values 
between. The curves may be regarded as the perpendicular 
projection of certain curves on the surface in space, the inter- 
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sections of the surface by equidistant planes parallel to th 
zy plane. 

The method is the generalization of the scale-representation 
of a function of one variable. For a relation between ¢ and x 
represented by a curve with ¢ as ordinate and x as abscissa, is 
transformed into a scale representation by perpendicularly 
projecting certain points of the curve onto the axis of 2x, the 
intersections of the curve by equidistant lines parallel to the axis 
of x and marking them with the value of ¢. A scale division in 
the case of a function of one variable corresponds to a curve in 
the case of a function of two independent variables. 

This method of representing a function of two independent 
variables by a plane drawing or we might also say of representing 
a surface in space by a plane drawing, is used by naval architects 
to render the form of a ship and by surveyors to render the form 
of the earth’s surface and by engineers generally. Let us apply 
the method to a problem of pure mathematics. 


The equation 
2+ pze+q=0 


defines z as a function of p and q. Let us represent this function 
by taking p and q as rectangular coérdinates and drawing the 
lines for equidistant values of z. 

For any constant value of z we have a linear equation between 
the variables p and q, and therefore it is represented by a straight 
line. This line intersects the parallels p = 1 and p= — 1 at 
the points g= — 22—z and g= —2#+2. Let us calculate 
these values for z = 0; + 0.1; + 0.2 --- + 1.3 and in this way 
draw the lines corresponding to these values of z as far as they 
lie in a square comprising the values p= —1 to +1 and 
qg=,—-1to +1. Fig. 44 shows the result. On this diagram 
we can at once read the roots of any equation of the third degree 


of the form 


2+ pe+q= 0, 
where p and q lie within the limits — 1to+1. Forp = 0.4 and 
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q = — 0.2, for instance, we read z = 0.37, interpolating the value 
of z according to the position of the point between the lines 
z= 0.3 and z = 0.4. We also see that there is only one real 
root, for there is only one straight line passing through the point. 
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Fic. 44, 


On the left side of the square there is a triangular-shaped region 
where the straight lines cross each other. To each point within 
this region corresponds an equation with three real roots. For 
example, at the point p= — 0.8 and ¢g= + 0.2 we read z= 
— 1.00; + 0.28; + 0.72. On the border of this region two roots 
coincide. 

For values of p and q beyond the limits — 1 to + 1 the diagram 


may also be used. We only have to introduce 2’ = z/m instead 


of 2 and to choose m sufficiently large. 
Instead of 


e+ oz+q=0 


: 
; 
7 
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we obtain 
msz!> + pmz’ + q=0, 
or dividing by m', : 
B Pp q 
a ar it ped gh 0, 
or 
3” =. p's! + q’ i 0, 
where 


P 
i reno ¢=4. 


By choosing a sufficiently large value of m, p’ and q’ can be 
made to lie within the limits — 1 to + 1 so that the roots 2’ 
may be read on the diagram. Multiplying them by m we 
obtain the roots z of the given equation. 

A function of two independent variables need not be expressed 
in an explicit form, but may be given in the form of an equa- 
tion between three variables 


g(u, v, w) = 0, 


and we may consider any two of them as independent and the 
third as a function of the two. The graphical representation 
may sometimes be greatly facilitated by modifying the method 
described before. The curves for constant values of one of the 
three variables, say w, are not plotted by taking w and v as 
rectangular codrdinates, but they are plotted after. introducing 
new variables x and y, x a function of wu and y a function of v and 
making x and y the rectangular codrdinates. 

In some cases, for instance, we can succeed by a right choice 
of the functions x = ¢(u) and y = y¥(v) in getting straight lines 
for the curves w= const. This will evidently be the case, 
when the equation g(u, 7, w) = 0 can be brought into the form 


a(w)e(u) + b(w)¥(v) + e(w) = 0, 


a, b, c being any functions of w, g any function of wu and y any 


function of 2. 
For introducing 
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t= (u), y= ¥() 
the equation will become 


ax+ by+c=0, 


where a, b, ¢ are constants for any constant value of w. 

As an example let us consider the relation between the true 
solar time, the height of the sun over the horizon, and the declina- 
tion of the sun for a place of given latitude. Instead of the 
declination of the sun we might also substitute the time of the 
year, as the time of the year is determined by the declination of 
the sun. Our object then is to make a diagram for a place of 
given latitude from which for any time of the year and any 
height of the sun the true solar time may be read. 

In the spherical triangle formed by 
the zenith Z, the north pole P (if we sup- 
pose the place to be on the northern 
hemisphere) and the sun S (Fig. 45), the 
sides are the complements of the decli- 
nation 6, the height A, and the latitude 
g. The angle t at the pole is the hour 
angle of the sun, which expressed in 
time gives true solar time. 

The equation between these four quantities may be written in 
the form 


Fia. 45. 


sin h = sin gsiné + cos ¢ cos 6 cost. 


The latitude ¢ is to be kept constant, so that #, h, 6 are the only 
variables. 
Now let us write 


x=cost, y=sinhk, 
so that the equation takes the form 
y = singsin6+ 2 cos g cos 6. 


When a and y are plotted as rectangular coérdinates, we obtain 


—— 
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a straight line for any value of 6. Let us draw horizontal lines 
for equidistant values of h = 0 to 90° and vertical lines for equi- 
distant values of ¢ = — 180° to + 180° or expressed in time 
from midnight to midnight (Fig. 46). In order to draw the 


x 


oe 
Spee 


° 
latitude — 41 


. 46. 


straight lines 6 = const., let us calculate where they intersect 
the vertical lines corresponding to «= —1 and x= +1 or 
expressed in time corresponding to midnight and to noon. For 
x = — 1 we have y = — cos (g + 5), and for 2 = + 1 we have 
y = cos (g — 6). Let us draw a scale on the vertical x = — 1 
showing the points y = — cos (y+ 4) for equidistant values of 
(g + 5) and a scale on the vertical x = + 1, showing the points 
y = cos (g — 6) for equidistant values of g — 6. The scale is 
the same as the scale for h, with the sole difference that the values 
of g — 6 are the complements of h and the values of gy + 6 the 
complements of — h. For a latitude of 41°, for instance, we 
have 


For 6 (i) seve y—sd 
RURITS Gi Meee siete eh ene ad esata aisha Goyal sie 23:0" 64.5° Leta 
September 23 and March 21........ 0 41° 41° 
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The values of ¢-+ 6 and yg — 6 furnish the intersections with 
the verticals 2 = — 1 and x = +1, so that the straight lines 
can be drawn corresponding to these days of the year. The two 
outward lines are parallel but the middle line is steeper. Their 
intersections with the horizontal line h = 0 show the time of 
sunrise and sunset. Strictly speaking the straight lines do 
not correspond to certain days. The straight line determined 
by any value of 6 changes its position continually as 6 changes 
continually. But the changes of 6 during one day are scarcely 
appreciable unless the drawing is on a larger scale. 


If in the equation 
ax-+ by+c=0 


a and 6b are independent of w, only ¢ being a function of w, all 
the straight lines w= const. are parallel. In this case we are 
not obliged to draw the 
straight lines w = const. 
It will suffice to draw a 
line perpendicular to the 
lines w= const. and a 
scale on it that marks the 
points corresponding to 
equidistant values of w. 
On the drawing we place a 
sheet of transparent paper 
Fic. 47. or celluloid,on which three 

straight lines are drawn is- 

suing from one point in the direction perpendicular to the u-scale, 
v-scale and w-scale (Fig. 47). If we move the transparent material 
without turning it and make the first two lines intersect the u-and-v 
scale at given points, the w-scale will be intersected at the point 
corresponding tc the value of w. This method has the advantage 


1 That is to say, the moment when the center of the sun would be seen on 
the horizon, if there were no atmospherical refraction. To take account of 


the refraction, the line h = — 0.6° would have to be considered instead of 
i 0) 
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that we can use the same paper for a great many relations of 
three variables, as we can place a great many scales side by side. 
Or, in the case of one relation only, we may divide the region of 
the values u,v, w into a number of smaller regions and draw three’ 
scales for each of them, placing all the u-scales or v-scales or 
w-scales side by side. The drawing will then have the same 
accuracy as a drawing of very much larger size in which there 
is only one scale for each of the three variables. 

§ 10. Depiction of One Plane on Another Plane.—Let us now 
-_ consider two quantities x and y each as a function of two other 
quantities w and v 


= g(u, v); 


y = o(u, 2). 


In order to givea geometrical meaning to this relation between 

- two pairs of quantities let us consider z and y as rectangular 
coérdinates of a point in a plane and uw, v as rectangular coérdi- 
nates of a point in another plane. We then have a corre- 
spondence between the two points. When the functions ¢(u, v) 

-and y¥(u, v) are defined for the values wu, v of a certain region, 
they will furnish for every point u,v of this region a point in 
the zy plane. Let us call this a depiction of the wv plane on 
the zy plane. Similarly a function of one variable x = ¢(u) 
might be said to depict the w line on the x line. We may there- 
fore say that the depiction of one plane on another plane is, in 
a certain way, the generalization of the idea of a function of one 
variable. Let us suppose ¢(u, v) and y(u, v) both to have only 
one value for given values of w and v for which they are defined. 
Then there will be only one point in the ay plane corresponding 
to a given point in the w plane. But to a given point in the 
ay plane there may very well correspond several points in the 
uv plane. 

Let us try to explain this by a graphical representation of the 
depiction of planes on each other. For this purpose we draw 
the curves x = const. and y = const. in the wv plane for equi- 
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distant values of 2 and y. In the zy plane they correspond to 
equidistant lines parallel to the axis of x and to the axis of y. 
The point of intersection of two lines = a and y = b corre- 
sponds to the points of intersection of the curves 


g(u,v) =a and ¥(u,v) = 8, 


in the wv plane. If in a certain region of the wo plane, that 
we consider, they intersect only once there is only one point in 
the region of the wv plane considered and one point in the zy 
plane corresponding to each other. Fig. 48 shows the depiction 
of part of the wo plane on part of the zy plane. We have a net 
of square-shaped meshes in the zy plane and corresponding is a 
net of curvilinear meshes in the wv plane. 

Let us consider the curves 2 = const. in the wv plane as the 
perpendicular projections of curves of equal height on a surface 
extended over that part of the wo plane. From any point P 
of the surface corresponding to the values u, » we proceed an 


0.2 0,3 0,4 0,5 0.6 


Fic. 48. 
infinitely small distance, w changing to wu + du, v to v+ dv and 
x to x+ dz, where 
oe we 
dz = ay ee + ae dv. 


Let us write 
du = cosads, dv = sin ads, 


where ds signifies the length of the infinitely small line from 
u, v to w+ du, v+ do in the wo plane and a the angle its direc- 
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tion forms with the positive axis of x. Let PN be a straight line 
whose projections on the w and v axis are equal to dy/du and 
d¢/dv and let us write 


r being the positive length of PN and } the angle between its 
direction and the positive axis of z. Then we have 


fe) 0 
i= qe t 5, = rds cos (a — d), 


or 


= r cos (a — A)... 


dx/ds measures the steepness of the ascent. It is positive when 
the direction leads upward and negative when it leads downward 
and its value is equal to the tangent of the angle of the ascent. 


From the equation 


a = + cos (a — di) 


we see that the ascent is steepest for a = \, where dz/ds = r. 
The line PN in the w, v-plane shows the perpendicular projection 
of the line of steepest ascent on the surface x = ¢(u, v) and the 
length of PN measured in the same unit of length in which w and 
” are measured is equal to the tangent of the angle of the ascent. 
Let us call the line PN the gradient of the function ¢(u, v) at the 
point u,v. The direction of the gradient is perpendicular to the 
curve ¢(u, v) = const. that passes through the point wu, v; for in 
the direction of the curve we have 


dz 
gee 


and therefore 
a—r= = 90°. 


If PN’ is the gradient of the function ¥(u, v) at the point wu, v, the 
angle between PN and PN’ must be equal to the angle formed. 
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by the curves « = const. and y = const. that intersect at the 
point u, v, or equal to its supplement according to the angle of 
intersection that we consider. 

Suppose the gradients PN and PN’ do not vanish in any of 
the points in the region of the wv plane that we consider and 
that their length and direction vary as continuous functions of 
wand v. Let us further suppose that the gradient PN’ (com- 
ponents: dy/du, dy/dv) is for the whole region on the left side 
of the gradient PN (components: d¢/du, d¢/dv), or else for the 
whole region on the right side of the gradient PN, then it fol- 
lows that any one of the curves x = const. and any one of the 
curves y = const. can only intersect once in the region considered. 

This may be shown by considering the directions of the curves 
x = const. and y = const. in the wo plane. Let us consider 
that direction on the curve y = const. in which z increases. If 
this direction deviates from PN the deviation must be less than 
90°, because dx/ds and therefore cos (a — ) is positive. Let us 
further consider that direction on the curve x = const. in which 
y increases. If it deviates from the direction of PN’ the devia- 
tion must be less than 90°. Let us call these directions the 
direction of 2 (on the curve y = const.) and the direction of y 
(on the curve 2 = const.). Now if the gradient PN’ is on the 
left of the gradient PN the y direction must also be on the left 
of PN (for if it were on the right of PN_ being perpendicular to 
PN it would form an obtuse angle with PN’) and therefore it 
must be on the left of the x direction (for if it were on the right, 
PN’ being perpendicular to the 2 direction would form an obtuse 
angle with the y direction, which we have seen to be impossible). 
Similarly it may be seen, that if PN’ is on the right of PN, the 
direction of y will also be on the right of the direction of x. If 
therefore PN’ is on the same side of PN in the whole region 
considered, the direction of y will also be on the same side of the 
direction of x for the whole region considered. This excludes 
the intersection of two curves x = const. and y = const. in more 
than one point. For, suppose there are two points of inter- 
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section and we pass along the curve y = const. in the direction of 
x. At the first point of intersection we pass over the curve 
x = const. from the side of smaller values of x to the side of 
larger values of x. Now if the values of 2 go on increasing 
as we go along the curve y = const. we evidently cannot get 
back to a curve 2 = const. corresponding to a smaller value of x. 
The only possibility of a second point of intersection would be 
that the direction in which the value of x increases on the curve 
y = const. becomes the opposite, so that in advancing in the 
same direction in which we came x would decrease again. 

The same holds for the curve 
x=const. If we pass from one 
point of intersection with a 
curve y = const. along a curve 
x = const. to a second point 
of intersection with the same 
curve the only possibility is 
that the direction of y also be- 
comes opposite. This is ex- 
cluded as in contradiction with Fic. 49, 
the direction of y being on the 
same side of the direction of a throughout the whole region (Fig.49) 

It will be useful to look at it from another point of view. Let 
us consider a point A in the wv plane corresponding to the 
values wu, v and let us increase wu and v by infinitely small positive 
amounts du and dv, so that we get four points ABCD, forming a 
rectangle corresponding to the coérdinates. 


A:u,v; B:ut+ dur; C:uvotd; D:u+ duvt+ d. 


In the zy plane these points are depicted in the points 4, 
B, C, D, the intersections of two curves u and u + du with two 
curves v and v» + dv (Fig. 50). 

The projections of the line AB in the zy plane on the axes of 
coordinates are obtained by calculating the changes of x and y 
for a constant value of v and a change du in the value of x 
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dx, = ss du, dy= oF dh 


Similarly the projections of AC are obtained by calculating the 
changes of x and y for a constant value of w and a change dv in the 
value of v 


Bik sels 
dx, = re dv, dip =a de, 


Denoting the lengths of AB and AC by ds; and ds. and the angles 
that the directions of 4B and AC form with the direction of the 


Fie. 50. 


positive axis of x (the angles counted in the usual way) by yj 
and y2 we have: 


dx, = ds, cosyi1, dy, = ds; sin y1 


and 

dx, = ds. cos y2, dy, = ds sin y, 
or 

0 ei Ok ee 

du Gat Cue ee 
and 

eer deh 0 an ee 

av edo on aa 
We may call 


e Bo dg 2 oy 2° 
du (35) + (*) 


the scale of depiction at A in the direction AB and 
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ds, _ N de\? , (avy? 

dv ~ Ge ay, 
the scale of depiction at A in the direction AC. It is here under- 
stood that the wv plane is the original, which is depicted on the 
zy plane. If we take it the other way the scales of depiction 
in the directions AB and AC are the reciprocal values du/ds, 
and dv/dsp. 

The area of the parallelogram ABCD in the zy plane is 


ds,dsz sin (y2 — 1) = (Gee _ sec) dud. 
According to the way in which the angles y. and y are defined 
sin (yz — 71) is positive, when the direction AC points to the left 
of the direction AB (assuming the positive axis of y to the left 
of the positive axis of x), and sin (yz — 71) is negative, when AC 
points to the right. Now dudv is equal to the area of the rectangle 
ABCD in the wo plane. Therefore the value of 


dg dy 00 Oy 
Ou Ov Ov OU 


is the ratio of the areas ABCD in the two planes and its positive 
or negative sign denotes the relative position of the directions 
AB and AC in the zy plane. We may call this ratio the scale 
of depiction of areas at the point A. 

dgdy dagdy 

au dv dv Ou 


is called the functional determinant of the functions ¢(u, v) and 
y(u, v). 

We have found the scale of depiction of lengths in the direc- 
tions AB and AC. Let us now try to find it in any direction 
whatever. From any point A in the wv plane, whose coérdinates 
are u and v, we pass to a point D close by whose coérdinates are 
u + Au, v+ Av. In the zy plane we find the corresponding 
points A and D with coédrdinates (Fig. 51). 
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= o(u, 2) _ a+ Ar = g(u+ Au, v+ Ar) 
At = yur) Phy + ay = Wut du, 0+ Ao) 


We expand according to Taylor’s theorem, and writing for 
shortness 


de a oe By ae ppg 
au £2 ant "8 aut 8. an 


Gu = 


we find 
Ax = ¢,Au + g,Av + terms of higher order, 


Ay = Au + W,Av + terms of higher order. 


: 
: 
: 
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The length of AD and the angle of its direction we denote by 
Ar and a in the wv plane and by As and X in the zy plane. 
The limit of the ratio As/Ar, to which it tends, when D approaches 
A without changing the direction AD is the scale of depiction 
at the point A in the direction AD. 
Writing 

Au = Ar cosa, 

Av = Ar sina, 
we obtain 


Ax = (gy CoS a+ ¢, sin a)Ar + terms of higher order, 
Ay = (Yu cos a + yy sin a)Ar + terms of higher order. 


Dividing by Ar and letting Ar decrease indefinitely, we have in 
the limit 


dx } 
pe Gu COSa+ gy Sina, 


» 
x eal 
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nm 


For de/dr and dy/dr we may also write ds/dr cos ), ds/dr sin }. 


ar °oS X= gucosa+ ¢g sina, 


ads. . 
Gp tn ® = Yucosa+ yy sina. 


These equations show the scale of depiction ds/dr corresponding 
to the different directions \ in the 2, y-plane and a in the u, v= 
plane. 

By introducing complex numbers we can show the connection 
still better. ; 2 
Let us denote " 


i: _ de dy. ds ~ 
Bat Be dee? 
a= gut yu, 
x = oy + Wl. 


Multiplying the second of the two equations by 7 and adding 
both they may be written as one equation in the complex form: 


Z2=2cosat wsina. 


The modulus of 2 is the scale of depiction of the wv plane at the 
point A in the direction a. The angle of z gives the direction in 
the zy plane corresponding to the direction a. For a = 0 we 
have z = 2 and for a = 90°, z = ~. 
Let us substitute 
i ett a e% . a Ct — emot 
Cove es gars li Sp 
and write 
z+ 2/4 zy — 2/t 
SS] b ve > aoe 
2 5: 2 


so that the expression for z becomes 
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z= ae“ + be, 


This suggests a simple geometrical construction of the complex 
numbers z for different values of a. The term ae“ is represented 
by the points of a circle described by turning the line that 
represents the complex number a round 
the origin through the angles a=0-- -2z. 
The term be~* is represented by the 
points of a circle described by turning 
the line that represents b round the ori- 
gin in the opposite direction through the 
angles a = 0--- — 27 (Fig. 52). The 
addition of the two complex numbers 
Fic. 52. ae“ and be™ for any value of @ is easily 
performed. The points corresponding 
to the complex numbers z describe an ellipse, whose two princi- 
pal axes bisect the angles between aand b. This is easily seen 
by writing 


a= ners b = melwtes, 


ao corresponds to the direction bisecting the angle between a 
and 6 and a; denotes half the angle between a and 6 (positive or 
negative according to the position of a and 6). 


s= rye (coat a)t + roe ora—a)é 
or 
eT aot << reir ent + Toe —(a—a)i 


(71 + 12) cos (a — a1) + (m1 — 72) sin (a — ay)t. 


Denoting the coérdinates of the complex number ze™™ by £ and 
we have 


E 


nee i: (a — a) and 


eager sin (a — a), 


and consequently the equation of an ellipse 
2 n° 
ai tp) a (ae 


Fe 
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This ellipse turned round the origin through an angle equal to 

ao gives us the points corresponding to z. The principal axes 

are 2(7. + ™) and 2(r;— m) (Fig. 53). The construction of 
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Fig. 53 is obvious. After plotting % and x we find ~/i and 
— x/t by turning AZ, through a right angle to the right and to 
the left. From these points lines are drawn to Z;._ The bisection 
of these lines give a and b. 

The figure shows that in case a and b have the same modulus, 
the triangle — Z,/i, Z,, Z2/i becomes equilateral and AZ; is per- 
pendicular to the line joining — Z2/i and Z2/t. In this case AZ, 
and AZ, would have the same or the opposite direction. But as 
Z = out Wul, % = or + Wot, this would mean that guy — gydu 
= 0. 

The radii of the ellipse (Fig. 53) measured in the unit used 
give the different scales of depiction corresponding to the dif- 
ferent directions in the xy plane. We might also say the ellipse 
is the image in the zy plane of an infinitely small circle in the 
uv plane, magnified in the proportion of the infinitely small radius 
to 1, with its center in A. 

Z, corresponds to a = 0 and Z; to a = 90° and for a = 0 to 90° 
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Z moves on the ellipse from Z, to Z, through the shorter way. 
— Z, corresponds to a = 180° and — Z, to a= 270°. Now we 
have shown above that a positive value of the functional deter- 
minant ¢,y%» — ¢v¥, means that Z, is on the positive side of Z1, 
so that in this case Z moves in the positive sense (that is, in the 
direction from the positive axis of z to the positive axis of y) with 
increasing values of a. With a negative value Z moves in the 
opposite direction. 

Let us now suppose that the curves 2 = const. and y = const. in 


the wv plane intersect except on a certain curve where their direc-. 


Fic. 54. 


tions coincide in the way shown in Fig. 54. On this curve the 
functional determinant D = gy, — gw, must vanish because 
the directions of the gradients coincide. Let us see what the 
depiction on the zy plane is like. 

Running along one of the curves y = const., say -y = 9, 
toward the curve D = 0 we intersect the curves 2 = a4, 23, 2 
until at the point A on the curve x = 2x, we reach the curve D = 0. 
In the xy plane the corresponding path is a parallel to the axis 
of x at a distance y passing through a4, 23, 2 and reaching a 
point A at 2. If we now proceed on the curve y = y; in the 
uv plane beyond the curve D = 0, we again intersect the curves 
%, 3, etc., but in the inverse order. Thus the corresponding 
path in the ay plane does not pass beyond A, but turns back 


. 


we a 
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through the same points 2, y1; 23, y1, etc. The same holds for 
any of the other lines y = const. If we trace the line in the 
zy plane that corresponds to the points in the wo plane, where 
the curves x = const. and y = const. touch, we find the depiction 
of the wv plane only on one side of the curve in the ay plane. 
The other side has no corresponding points u, v. However to 
every point C on this side of the curve, there are two correspond- 
ing points C in the wv plane, one on either side of the curve 
D=0. Imagine two sheets of paper laid on the ay plane; let 
them both be cut along the curve AB. Retain only the two 
pieces on this side of the curve and paste them together along 
the curve. The wo plane is in this way depicted on the paper 
in such a way that there is one point and one only on the paper 
corresponding to each point in 
the region of the wo plane con- 
sidered. The curve D=0 in 
the wv plane corresponds to the 
rim where the two pieces of pa- 
per are pasted together. Any 
line straight or curved passing 
over the curve D = 0 in the wo 
plane,corresponds to a line running from one of the sheets onto the 
other. It need not change its direction abruptly when it reaches 
the rim and passes onto the other sheet. For it may touch the 
rim in the direction of its tangent. This is actually the rule 


and the abrupt change of direction is the exception. Any line 


LAL (Fig. 55) in the wv plane, whose tangent as it crosses the 
curve D = 0 at A does not coincide with the common tangent 
of the curves x = const. and y = const. will correspond to a line 
in the zy plane, that does not change its direction abruptly 
when it touches the rim. 

This is best understood analytically. Let us consider corre- 
sponding directions at the points A in the wo plane and in the 
ay plane. We have seen above that corresponding directions 
(Fig. 56) are connected by the equations 
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dz ; 
cosAT HS T= cos a sin 
dr dec + % ae 


d 
sinh 7. = 7 Yu cosa+ yp, sina. 


At the point A we have 
Puby — ovbu = 0. 


Assuming that the gradients at A do not vanish, so that we 
can write 


Gu=T COSY, @ =rsiny, 


Vu 


where r and 7’ are positive quantities, the equation ¢u»— g»t,=0 
reduces to sin (y — y’) = 0, that is, y= y’ or y = y’ + 180°. 
It follows therefore that: 


r cosy’, =r siny’, 


es 
cosh 7 = 7 cos (a—y¥), 


by ds , , 
sin Aa = r’ cos (a — 7’) = =r’ cos (a — y). 
Consequently for all directions a in the wv plane for which 
cos (a — 7) is not zero, we have 
r! 
tg A= = Ho 


ee ae 
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That is to say, we have in the zy plane only one fixed direction 


\ and the opposite corresponding to all the different directions 


a except only a direction for which cos (a — y) = 0. In the 
latter case, that is, when the direction @ is perpendicular to the 
direction y of the gradient, 7. e., in the direction of the curves 
x = const. and y = const., we have 


ds 
cos ar = 0, 


ae 
sin Ae 0. 
Therefore ds/dr = 0 and \ remains indeterminate. Any direction 
d for which tg \ differs from + r’/r corresponds to a fixed direction 
a= 7+ 90° ora = y — 90°, while ds/dr = 0. 

As the curve D = 0 is depicted on the rim of the two sheets 
of paper, all those lines that intersect the curve D = 0 in a 
direction different from the direction of the curves x = const. 
and y = const. are depicted in the zy plane as curves having 
their tangent at A in common with the rim. All lines in one of 
the sheets of paper that touch the rim at A in a direction differ- 
ent from that of the rim must be the depiction of lines in the wv 
plane that reach A in the direction of the lines x = const. and 

= const. The scale of depiction is zero in the direction of the 
curves xz = const. and y= const. In any other direction a 
we find it different from zero for: 


.. = V(r +r”) cos’ (a — 4). 


It is a maximum in the direction a = y or y + 180° perpendicular 
to the curves x = const. and y = const. 

It may help to understand all these details if we discuss an 
example where the depiction of the wo» plane on the zy plane 
has a simple geometrical meaning, the planes being ground plan 
and elevation of a curved surface in space. The rim in the 
ay plane is the outline of the surface, the projection of those 


80 GRAPHICAL METHODS. 


points where the tangential plane is perpendicular to the plane 
of elevation. 

Suppose a cylinder of circular section cut in two half cylinders 
by a plane through its axis. Suppose one of the half cylinders 
in such a position that its axis 
forms an angle 6 with the 
ground plan, the plan of ele- 
vation being parallel to its — 
axis, Fig. 57. Let us intro- 
duce rectangular coédrdinates 
u,v in the ground plan and 
rectangular codrdinates 2, y 
in the plan of elevation. A 
point P on the cylinder is de- 
fined by certain values u, v 
which define its ground plan 
and certain values x, y which 
define its elevation. It is 


a BT. easily seen from Fig. 57 that 
we have 
-£=4 
and 
= fae Lien 
Rapist gS va 


where a is the radius of the section. Now let us consider the 
elevation of the points P as a depiction of their ground plan. 
The functions ¢(u, v) and y(u, v) in this case are 
elu, 0) = u, | 
1 Mt 
u, 0”) =ut —Ve@— 
Yu, v7) = 4 gon eke v, 
and 
v 


“ee ee 
cos 6Va? — 2 


Cia il Py = 0; Vu = tg 6, Wy =— 


v 
cos 6V. a2 —: x 


Pub» — Pu oe 
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The functional determinant vanishes for » = 0 on the line EF. 
The lines y = const. are the intersections of the cylinder with 
horizontal planes. In the plan of 
elevation they are straight hori- 
zontal lines; in the ground plan 
they are ellipses (Fig. 58). As we 
pass along one of these curves we 
cross the line EF in the ground 
plan but we only touch it in the 
_ plan of elevation; retracing the hori- 
zontal line back again. The lines 
x= const. are straight lines in both 
planes, but in space they corre- 
spond to ellipses. Again as we 
cross EF in the ground plan we 
only touch it in the plan of eleva- | 
tion and retrace the vertical line down again. Any curve on 
the cylinder that crosses EF in a direction not perpendicular to 
the plan of elevation is projected in the plan of elevation with 
EF as its tangent. For the real tangent in space lying in the 
tangential plane of the cylinder can have no other projection, if 
not perpendicular to the plan of elevation. In this latter case 
the projection of the tangent is a point 
and the tangent of the elevation is deter- 
mined by the inclination of the osculatory 
plane. 

There is a particular case to be consid- 
ered, when the curve D = 0 in the wo plane 
coincides with one of the curves 2 = const. 

Fic. 59. or y = const. (Fig. 59), assuming the gra- 

dients of the functions g(u, v) and (uw, v) 

not to vanish at the points of this curve. We have seen that at 

a point where D = 0 the scale of depiction must vanish in the 

directions of the curve z = const. or y= const. Let the curve 

D = 0 coincide with a line x = const., then it follows that the 
7 


Fig. 58. 


D=o0 


\ 
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length of the depiction of this curve is zero and the depiction 
must be contracted in a point. For the length of the depiction 
of a curve x= const. 1s given by an integral 


ds 
{Ga 


where dr denotes an element of the curve and ds/dr the scale 
of depiction in the direction of the curve. As ds/dr is zero all 
along the curve the integral must necessarily vanish. 

As an example let us con- 
sider 


The lines x = const. in the wv 
plane are equilateral hyper- 
bolas, the lines y = const. are 
parallels to the axis of wu (Fig. 
60). Along the axis of uw we 
eed l pees have at the same time y = 0, 
et ene x=0 and D=v=0. The 
whole axis of w is depicted in 

the point x = 0, y = 0 of the zy plane. 

Let us finally consider the case where the scale of depiction 
at any point is the same in all directions, though it need not be 
the same at different points. 

Writing as before 


ala Pu Ss Wut, 2 Yv + Vol, 
ds 


_ dz dy. - 
i ipa dea ile 


the connection between the scale of depiction ds/dr and the 
angles \, w determining corresponding directions in the zy plane 
and in the wv plane is given by the equation 


2 = 2cosa+ xsina, 


Bae bi 
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or 
a= ae* + be 
where 
a=}ata/i), b= 2G — aii). 


In the case where the scale of depiction ds/dr, that is to say, the 
modulus of z, is independent of a, one of the constants a or b 


must vanish, as we see at once from the construction of z (Fig. 
52). Let us consider the case b = 0, 


bee he 
let a HS eM, 
The complex number a may be written |a| e, where | a | 
denotes the modulus of a and ap the angle. Both may vary 
from point to point, but at every point they have fixed values. 
Consequently we have 
—_ la| and AX\=atay 
That is to say, from an angle a determining a direction in the 
uv plane, we find the angle determining the corresponding 
direction in the zy plane by the addition of a fixed value apo. 
Any two directions a, a’ will therefore form the same angle as. 
the corresponding directions \, \’ in the zy plane. The same is 
true when a = 0 and z = be. The only difference is that in 
this latter case the direction of z rotates in the opposite sense 
with increasing values of a. 
Analytically depictions of this kind are represented by func- 
tions of complex numbers, 


atyi=fiutv) or «a+y=f(u— nm). 


Assuming the function to possess a differential coefficient we have 


Cts Cy, ‘ ; 
= — oe —4 = = 
edie aut Te 7); 


on ey. ; <3 
<=) —— —4-= = = y 
re + ant fu + n)t, 
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and therefore either 
a=2i or 4= — &/t. 
Hence in the first case 
a=iiaatali=u, b= }(a—a/i) =0 


and in the second case 
a=0, b=. 

§ 11. Other Methods of Representing Relations between Three 
Variables.—The depiction of one plane on another may be used 
to generalize the graphical representation of a function of two 
variables or a relation between three variables, as we prefer 
to say. 

As we have seen before, an equation 


g(x, y, 2) = 0 


between three variables 2, y, z can be represented by taking x 
and y as rectangular coérdinates and plotting the curves z = 
const. (Fig. 61) for equidistant val- 
ues of z. Suppose now the zy plane 
to be depicted on another plane. 
The lines x = const.,y = const. and 
z = const. will be represented by 
three sets of curves. The fact that 
a three values z, y, 2 satisfy the equa- 
tion g(x, y, 2) = 0 is shown geo- 
metrically by the intersection of 
the three corresponding curves in 
one point. , 


y=const. 


Another method for representing certain relations between 


three variables wu, v, w consists in drawing three curves, each 
curve carrying a scale. The values of u, v, w are read each on 
one of the three scales. The relation between three values w, 2, 
w is represented geometrically by the condition that the corre- 
sponding points lie on a straight line (Fig. 62). This method is 


Pe | oe 
’ r 
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far more convenient than the one using three sets of curves. It is 
less trouble to place a ruler over two points u, v of two curves 
and read the value w on the scale of the third than to find the 
intersection of two curves w = const. and v = const. among sets 
of others, pick out the curve w = const. that passes through the 


Fia. 62. 


same point and read the value of w corresponding to it. For we 
must consider that the curves corresponding to certain values 
of w and v are generally not drawn, but must be interpolated and 
so must the curve w = const. It is true that interpolations are 
necessary with both methods, but the interpolation on scales 
like those in Fig. 62 is easily done. 

It must however be understood that while the three sets of 
curves form a perfectly general method for representing any rela- 
tion between three variables, the other method is restricted to cer- 
tain cases. In order to investigate this subject more fully we 
shall have to explain the use of line coérdinates. 

When we apply rectangular coérdinates 2, y to define a certain 
point in a plane, we may say that 2 determines one of a set of 
straight lines (parallel to the axis of ordinates) and y determines 
one of another set of straight lines (parallel to the axis of abscissas) 
and the point is the intersection of the two (Fig. 63, J). A 
similar method may be used to determine a certain straight line 
in a plane. Let 2 determine a point on a certain straight line, 
x being its distance from a fixed point A on the line measured 
in a certain unit and counted positive on one side and negative 
on the other. Let y define a point on another straight line 
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parallel to the first, y being its distance from a fixed point B on 
the line measured in the same way as xz. The straight line 
passing through the two points is thus determined by the values 
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x and y and for all possible values of 2 and y we obtain all the 
straight lines of the plane except those parallel to the lines on 
which 2 and y are measured. For simplicity we choose AB 
perpendicular to the two lines (Fig. 63, IZ). Let us call x and y 
the line codrdinates of the line connecting the two’points x and 
y in Fig. 63, IJ, in the same way as 2 and y in Fig. 63, J, are 
called the point coérdinates of the point where the two lines 
x and y intersect. 
A linear equation between point codrdinates 


y= met pu 


is the equation of a straight line. That is to say, all the points 
whose coédrdinates satisfy the equation lie on a certain straight 
line. If, on the other hand, we regard 2 and y as line coérdinates 
we find the analogous theorem: all the straight lines whose 
line coérdinates satisfy the equation 


y=me+p 


pass through a certain point. The equation is therefore called 
the equation of the point. 


In order to show this let us first draw the line = 0, y = 
(APO in Fig. 64). If now for any value of we make AR = x 


aan 3 OT 
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and PQ = mz, the point of intersection of RQ and AP must be 
independent of x, for 


The ratio PO/AO determines the position of 0 and as it is 
independent of x and the positions of A and P are also inde- 
pendent of xz, the same is true for 0. 
For negative values of m, PO and 
AO have opposite directions so that 
O lies between A and P. 

For a given point O, we can find 
the corresponding values of mand u 
by joining O with the points A and 
the point corresponding to z= 1. 
If P and Q are the intersections of 
these lines with the line on which y 
is measured, we have BP = » and PQ= m. Any point in the 
plane thus leads to an equation 


y= m+ ph, 
except the points on the line on which 2 is measured. For 
m = 0 the equation reduces to 
Y=, 

that is, the equation of a point on the line on which y is measured. 

Instead of y = mz+ yu, we might also write x = m'y+ w’, 
and go through similar considerations changing the parts of x 
and y. This form does not include the points on the line on 
which y is measured, but it does include the points on the line 
on which 2 is measured. For these we have m’ = 0. 

The general equation of a point in line co6rdinates is given in 


the form 
ax + by+c= 0, 


from which we may derive either of the first-mentioned forms 
dividing it by a or b. 
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Dividing by ¢ another convenient form is obtained, 


ax by 


or writing 


x y 

past ARE te cces, 

i 5 e , 
ao determining the point of intersection of the line BO (Fig. 64) 
and the a-line, while y) determines the point of intersection of 
the line AO with the y-line. 


A curve may be given by an equation 
a(u)e + bi(u)y + a(u) = 0, 


in which a;(w), 6:(w), c:(u) are functions of a variable wu. Any 
value of w furnishes the equation of a certain point and as wu 
changes the point describes the curve. Let us suppose the curve 
drawn and a scale marked on it giving the values of w in certain 
intervals sufficiently close to interpolate the values of u be- 
tween them. ‘Two other curves are in the same way given by 
the equations 


a(v)e + (vy + @(v) = 0, 
as(w)a + b3(w)y + c3(w) = 0, 


and scales on these curves mark the values of v and w. 

Now we are enabled to formulate the condition which must be 
satisfied by the values w, v, w in order that the three corresponding 
points lie in one straight line. If 2 and y are the line codrdinates 
of the line passing through the three points, x and y must satisfy 
all three equations simultaneously. 


Consequently the determinant of the three equations must 
vanish 


ay (boc3 = bsc2) = aa(b3cr re by¢3) - a3(bic2 a bc) = 0, 


and, vice versa, if the equation between wu, 7, w may be brought 
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into this form where a, by, ¢: are any functions of wu, a, be, ¢ any 
functions of v and as, 63, cs any functions of w, we can form the 


equations 
nx + by +c = 0, 


age + boy + & = 0, 
asx + bsy + cs = 0, 


and represent them graphically by curves carrying scales for 

u, v, w. The relation between u, 2, w is then equivalent to the 

condition that the corresponding points on the three curves lie 

on a straight line. But it must be remembered that only a 

restricted class of relations:can be brought into the required form, 

so that the method cannot be applied to any given relation. 
The equation of a point 


ax+ by tc=0 


remains of the same form, when the units of length are changed 
forwand y. If 2’ denotes the number measuring the same length 
as the number z but in another unit, the two numbers must have a 
constant ratio equal to the inverse ratio of the two units. There- 
fore, by changing the units independently, we have 


2=)r’, y= py’, 
and the equation of the point may be written 
az’ + buy’ + ¢ = 0, 
a’z’ + b'y’ +c =0, 


where a’ = da and J’ = wb. 

It is sometimes convenient to define the line coérdinates in 
another way. Let & and 7 denote rectangular codrdinates 
measured in the same unit, then the equation of a straight line 
can be written 


or 


n = tg v& + no, 


where ¢ is the angle between the line and the axis of £ and 1p, 
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the ordinate of the point of intersection with the axis of ». 
Now let us call tg g and 1 the line coérdinates of the straight 
line represented by the equation and let us denote them by x 
and y. Thus the values of x and y define a certain straight line 
and any straight line not parallel to the axis of ordinates may 
be defined in this manner. The condition that a straight line 
2, y passes through a point £, 7 is expressed by the equation 


n=wét+y, 
or 
y = —4e +4. 


If we fix the values of x and y, all the values &, 7 that satisfy this 
equation represent the points of the straight line z, y and we 
therefore call it the equation of the straight line. If, on the 
other hand, we fix the values of £ and 7, all the values z, y that 
satisfy the equation represent the straight lines that pass through 


the given point é, 7, and therefore we call it the equation of the ~ 


point. 
The more general form 
ax+ by +c=0 
can be reduced to 
a c 
a 


It therefore represents the equation of the point, whose rec- 
tangular codrdinates are £ = a/b and » = —c/b. The case 
where b = 0 or 

ax+e=0 


represents the equation of a point infinitely far away in the 
direction ¢ or the opposite direction ¢ + 180°, ¢ being defined by 


c 
Ee des ee ox 


All the straight lines, whose coérdinates 2, y satisfy the equation 


ax+e=0 


Ee 


~ 
; 
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correspond to the same value of x but to any value of y.- That 
is to say, they are all parallel and all the straight lines of this 
direction belong to them. 

Let us now discuss some of the applications of line codrdinates 


to the graphical representation of relations between three 
variables. 


The relation 


Be 


o fF OA aI MSO 


w= Ww 


eS 


may be written in the form 


log u + log v = log w, 
or 
x+y = log w, 
when 


zx=logu and y= log». 


Let us plot z and y as line co- 
érdinates on two parallel lines (Fig. 
65), with scales for the values of u 
and v. The equations x = log u 
and y = log v may be regarded as the equations of the points of 
these two scales. The equation 


Fig. 65. 


a+ y = logw 


for any value of w is the equation of a point. It can easily be 
constructed as the intersection of any lines 2, y satisfying its 
equation. For instance, the line x = log w, y = 0 and the line 
x= 0, y=logw. The first line is found by connecting the 
scale division uw = w of the u-scale with the point B, the second 
by connecting the scale division » = w of the v-scale with the 
point A. If the units of 2 and y are taken of the same length, the 
point of intersection will lie in the middle between the two lines 
carrying the u and » scales on a line parallel to the two other lines 
and the w-scale will be half the size of the other two (Fig. 65). 
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The relation 
Uv = Ww 
or 
log u + log v = log w 
expresses the condition that the three equations 
zx=logu, y=logr, «+y=logw 


are satisfied simultaneously by the same values of x and y, that 
is to say, that the three points on the w, v, w scales corresponding 
to the values of wu, v, w lie on the same straight line 2, y. 

The more general relation 


uv? = w, 
where a and @ are any given values, can be treated in the same 
manner. Thus the pressure and volume of a gas undergoing 
adiabatic changes may be represented. In this case we have 
pu = w, 
where p denotes the pressure, » the volume and k and w con- 
stants. 
For a given gas k has a given value, but w depends on the 
quantity of the gas considered. 
We write 
x=logp, y= log». 


The relation then takes the form 


x + ky = log w, 


and represents a point which may be con- 
structed by the intersection of any two 
straight lines x, y, whose codrdinates sat- 
isfy the equation, for instance 


x=logw, y=0 
and 


1 
x= 0, y = 7, log w. 


; 
1 
| 


———— 


GRAPHICAL REPRESENTATION OF FUNCTIONS. 93 


The first line connects the point B (Fig. 66) with the scale 


division p = w of the p scale and the second line connects the 
point 4 with the scale division of the v scale for which y = hk log w. 
A perpendicular from the point of intersection on AB meets it in 
O’ and as the ratio A0’/O’B is 
equal to the ratio of the seg- 
ments on the p and 2 scales 
log w/k log w = 1/k it is inde- 
pendent of w. All the points 
corresponding to different val- 
ues of w lie on the same par- 
allelto the p and v scales and 
the w scale may be obtained 
by a central projection of the 
p scale on this parallel from 
the center B (Fig. 66). We 


+10 


might dispense with the con- —1.0 
struction of the w scale as 
long as the straight line for 

: f ~1.5 
the w scale is drawn. For in 
using the diagram we gener- 
ally start with values po, vo ~2.0 


and want to find other values 
p, %, for which 


pk = Doro". Fig. 67. 


The straight line connecting the scale divisions p and »v intersects 
the w scale at the same point as the straight line connecting the 
scale divisions po and %, so that we need not know the value of 
povo®. It suffices to mark the point of intersection in order to 
find the value of p, when v is given or the value of v when p is 
given. 

Another example is furnished by the equation 


wtawty=0. 
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If we regard a and y as line codrdinates any value of w determines 
the equation of a point. We plot the curve formed by these 
points with a scale on it indicating the corresponding values of w. 
Any values of 2 and y determine a straight line whose inter- 
sections with the w scale furnish the roots of the equation. Each 
point of the w scale may be constructed by the intersection of 
two straight lines, whose coérdinates 2, y satisfy the equation, 
for instance 
z=0, y=—w and t=—w, y=0/ 


In Fig. 67 the w scale is shown for the positive values w = 0 to 
w= 2.5. 
In the same manner a diagram for the solution of the cubic 
equation 
w+awty=0. 
or of any equation of the form 


w+ae*+y=0 
may be constructed. 


$12. Relations between Four Variables —The method can be 
generalized for relations between four variables. 
Suppose four variables w, v, w, t are connected by the equation 


g(u, v; W,; t) ma 0, 


and let us assume that for any particular value ¢ = fy the resulting 
relation between w, v, w can be given by a diagram of the form 
considered consisting of three curves carrying scales for wu, » and 
w. Let us further suppose that for other values of ¢ the scales 
for w and v remain the same, but the scale-for w changes. Then 
we shall have a set of w scales corresponding to different values 
of t. Connecting the points that correspond to the same value 
of w we obtain a network of curves ¢t = const. and w = const. 
(Fig. 68). Any two values w, v furnish a straight line intersecting 


1 For small values of w, this combination is not good because the angle of 
intersection is small. One might substitute = 2, y = — w® — 2w for the 
first line. 


ore oe, Tt asi tnenincanteiaiaiiladetaliall 
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the network of curves. The points of intersection correspond to 
values of ¢ and w that satisfy the given relation. 
Any relation of the form 


may be represented in this way, g(u) denoting any function of 
u, Y(v) any function of v and f(t, 
w), g(t, w), h(t, w) any functions 
of ¢t and w. 

In this case we need only in- 
troduce the line coédrdinates 2, 
y, writing 


t= 9(u), y= (0). 


We then obtain a linear equation 
between z and y, 
“001-2 8 4B SE 


ft, wx + gt, w)y + hG, w) = 0, 


which for any given values of t ™ Wis 6. 

and w represents the equation of 

apoint. For a given value of ¢ and variable values of w we obtain 
a curve t = const. carrying a scale for w and for a series of values 
of ¢ we obtain a set of curves t = const. Similarly for a given 


- value of w and variable values of ¢ the equation furnishes a curve 


w = const., carrying a scale for ¢ and a series of values of w 
furnishes a set of curves w = const. From any given values 
of w and v the line coérdinates z and y are calculated and the 
points where this straight line defined by a and y intersects 
the netwerk of the curves ¢ = const. and w = const. furnish 
the values ¢, w that satisfy the relation together with the given 
values of w and v. The relation between the height, azimuth, 
declination of a celestial body and the latitude of the point of 
observation may serve as an example. Let h, a, 6 denote the 
height, azimuth and declination and ¢ the latitude. The angles 
1/2 — o, 7/2—h, 7/2 —6 are the three sides of a spherical 
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triangle PZS (Fig. 69) formed by the pole P, the zenith Z and the 
celestial body S. The azimuth is defined as the supplement of 
the angle PZS. 

The equation is 

sin 6 = sin g sinh — cos g cosh cosa. 


We write 


Fic. 69. 


z=cosa, y= sind, 
so that the equation becomes 
y = sin gsinh — xcos g cosh. 


We shall in this case use the second system of line codrdinates 
where x is the slope of the line measured by the tangent of the 
angle formed with the axis of 
abscissas and y is the ordinate 
of the intersection with the axis 
of ordinates. If £, » denote the 
rectangular codrdinates of the 
point, the equation of the points 
takes the form 


Mace Or =p er, 
so that in our case we have 
~=cos¢gcosh, n=singsinhk. 


The curves ¢ = const. and h = 
const. can be drawn by means 
of these formulas. It is easily 
seen that they are ellipses and 
that the curves y = const. are 
the same as the curves h = const. 
For a definite value of g and a 
variable value of h we find 


Fig. 70. 


cos” sin? 9 x 


) 


¢ = *"s Wes 
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: and for a definite value of 2 and a variable value of ¢ 


- Sy Ue 
cosh sin? hh = 
Any of the ellipses intersects all the others and in this way they 
form a network. A point of intersection of the ellipse g = ¢ 
and the ellipse h = ¢ also corresponds to the values h = c and 
¢ = %, as the ellipse g = ¢ is identical with the ellipse h = ¢ 
and g = & identical with h = c (Fig. 70). The easiest way to 
find this network consists in drawing the straight lines 


E+ 7 = cos (g — h), 
and perpendicular to them the straight lines 


—~— n= cos(¢+ fh), 


for equidistant values of g +h and g—h. The ellipses run 
diagonally through the rectangular meshes formed by the two 
systems of straight lines. The scales for g and h are written 
on the axis of codrdinates, both scales being available for both 
variables. The scale for 6 is written on the axis of ordinates 
and is identical with the scale for ¢ and h on this axis. For the 
ordinate corresponding to a given value 6 = cis sinc, and this is 
also the ordinate of the point where the ellipse g=c or h=c 
intersects the axis of ordinates. The scale for the azimuth cannot 
be laid down in exactly the same way as that for yg, h and 6 
because cosa determines the slope of the straight line a, y. 
Let us draw a parallel to the axis of ordinates through the point 
= 1, 7 = 0 and mark a scale for the azimuth on it, making 
n = cosa (Fig. 70). A line connecting the origin with any scale 
division of this scale has the slope of the line z = cos a, y = sin 6. 
To bring it into the position of the line 2, y it must be moved 
parallel to itself, until its point of intersection with the axis of 
ordinates coincides with the scale division 6. This suggests 
another way of using the diagram. Let a pencil of rays be 
drawn from the origin to the scale divisions of the azimuth scale 
(Fig. 70), and let it be drawn on a sheet of transparent paper 
8 
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placed over the drawing of the ellipses. For any given value 
of 6 it is moved up or down as the case may be so that the center 
of the pencil coincides with the scale division 6. As long as the 
celestial body does not materially alter its declination the dia- 
gram in this position will enable us to find any of the three 
values , h, a from the other two. 

As a second example let us consider the relation between the 
declination 6, the azimuth a, the hour angle ¢ of a celestial body 
and the latitude ¢ of the point of observation. 

The relation is found by eliminating the height f from the 
equation 

sin 6 = sin g sin h — cos g cosh cosa. 


For this purpose we express sin h and cos h by the other angles 


and substitute these expressions for sin h and cos h. 
We have 
cos h = cos 6 sin t/sin a, 


sin h = sin ¢siné+ cos ¢ cos 6 cost. 


Substituting these values we find 
sin 6=sin? ¢ sin 6+sin ¢ cos ¢ cos 6 cos t—cos ¢ cos 6 sin t ctg a, or 


cos” ¢ sin 6 = sin ¢ cos ¢ cos 6 cos t — cos ¢ cos 6 sin t ctg a. 


Dividing by cos? ¢ cos 6 we finally obtain 


sin ¢ 
tg 6 = tg gcost — : 


ctg a. 
COs ¢ 8 


In order to represent this relation graphically we introduce line 
coordinates 
x=ctga and y= tgéd 
and find 
sin ¢ 
=t i 
y AE y CON aes 


Let us use the second system of line codrdinates. The rec- 


tangular coérdinates £, 7 of the point represented by the equation 
are found from it equal to: 


ver) ll i 


~ the negative part of the axis and ¢ is 
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sin t 
é = 
cos ~ 


» n=tg gcost. 


The curves ¢ = const. are ellipses, 


The curves t = const. are hyperbolas, 


2 n° 
sin?t  cos?t 


The ellipses and hyperbolas are confocal, the foci coinciding 
with the points = + 1, 7 = 0, so that the curves intersect at 
right angles. 

The scale for g may be written on the axis of ordinates at the 
points where it intersects the ellipses. It is identical with the 
scale for 6, the ordinate in both cases 
being the tangent of the angle with the 
only difference that 6 is negative on 


not. The scale for t may be written 
on one of the ellipses corresponding to 
the largest value of ¢ that is to be taken 
account of. ‘This ellipse forms the 
boundary of the diagram, so that 
larger values of g are not. represented. 
Corresponding to the azimuth we draw 
a pencil of rays on a sheet of trans- 
parent paper, which is laid on the draw- 
ing of the curves. The center of the 
pencil is placed on the scale division 6 
and the azimuth is equal to the angles 
that the rays form with the positive direction of the axis of or- 
dinates (Fig. 71). It suffices to draw the curves and the rays 
only on one side of the axis of ordinates. At the apex of the 


a: pencil aie rays is is in Hileg thea rays | fee 
F ap “Tines of the hyperbolas, a = 15° comesponding to, = 
tot = 2" and so on. © 
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CHAPTER III. 


Tue GrapHicaAL Metuops oF THE DIFFERENTIAL AND 
INTEGRAL CALCULUS. 


§ 13. Graphical Integration—We have shown how the ele- 
mentary mathematical operations of adding, subtracting, multi- 
plying and dividing and the inverse operation of finding the 
root of an equation can be carried out by graphical methods and 
how functions of one or more variables may be represented and 
handled. But the graphical methods would lack generality and 
would be of very limited use, if they were not applicable to the 
infinitesimal operations of differentiation and integration. In- 
deed it is here that they are found of the greatest value. In 
many cases, where the calculus is applied to problems of natural 
science or of engineering, the functions concerned are given in a 
graphical form. Their true analytical structure is not known 
and as a rule an approximation by analytical expressions is not 
easily calculated nor easily handled. In these cases it is of vital 
importance that the operations of the calculus can be performed, 
although the functions are only given graphically. 

Let us begin with integration, because it is easier than differ- 
entiation and of more general application. 

Suppose a function y = f(x) given by a curve whose ordinate is 
y and whose abscissa isa. The problem is to find a curve, whose 
ordinate Y is an integral of the function f(z), 


Y= il Winds. 


Let us assume the unit of length for the abscissas independent 
of the unit of length for the ordinates. The value of Y measures 
the area between the ordinates corresponding to a and 2, the 
curve y = f(z) and the axis of x in units equal to the rectangle 


formed by the units of x and y. 
101 
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In the simple case where f(x) is a constant the equation 
y = f(x) = cis represented by a line parallel to the axis of a and 


y={ edz = c(x — a). 


Y is the ordinate of a straight line intersecting the axis of x at 
the point x =a. The constant ¢ is the change of Y for an 
increase of 2 equal to 1. 
If P is the point on the 
axis of x fora = — 1 and 
Q the point where the line 
y = c intersects the axis 
of ordinates (Fig. 72) the 
desired line is parallel to 
Fic. 72. PQ. It is constructed by 

drawing a parallel to PQ 

through the point 2 = a on the axis of x (Fig. 72, where a = 0). 
When a given value ¢; is added, so that the equation becomes 


Y=c(a—a)+aQ 
it amounts to the same as when the straight line is moved in the 


direction of the axis of ordinates through a distance ¢. For 
x = a we then have Y = «q, so that we obtain the line 


Y=c(e—a)+ 4, 
by drawing a parallel to PQ through the point 2 = a, y = q. 
In the second place let us assume that the line y = f(x) consists 
of a number of steps, that is to say, that the function has different 


constant values in a number of intervals x = 21 to a, a to 23, 


etc., while it changes its value abruptly at a, 23, etc. The 
line presenting the integral 


Y= ip f(x)dx 


does not change its ordinate abruptly. It consists of a con- 
tinuous broken line, whose corners have the abscissas Xe, X3, etc. 


eee ia 
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The directions of the different parts are found in the way just 
described by the pencil of rays from P to the points a, B, ¥, ete. 
(Fig. 73), where the horizontal lines intersect the axis of ordinates. 


Fig. 73. 


To construct the broken line we draw a parallel to Pa through 
the point z = 2 (in Fig. 73 2 is equal to 0) as far as the vertical 
2= 2. Through the point of intersection with the vertical 
x = a we draw a parallel to P8 as far as the vertical x = ay. 
Through the point of intersection with the vertical x = x3 we 
draw a parallel to Py and so on. 

Finally let us consider the case of an arbitrary function y = f(z) 
represented by any curve. In order to find the curve 


Y= f fede 


we substitute for y = f(x) a function consisting of different 
constant values in different intervals and changing its value 
abruptly when a passes from one interval to the next, so that 
the line representing this function consists of a number of steps 
leading up or down according to the increase or decrease of f(z). 
These steps are arranged in the following way. The horizontal 
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part A,A» of the first step (Fig. 73) starts from any point Ai 
of the given curve. The vertical part 4,B, and the following 
horizontal part B,B, are then drawn in such a manner that BB, 
intersects the curve and that the integral of the given function 
as far as the point of intersection K, is equal to the integral of the 
stepping line as far as the same point. That is to say, the areas 
between the stepping line and the given curve on both sides of 
the vertical part 4.B, have to be equal. When K; is fixed the 
right position of A,B; may be found by eye estimate. The eye 
is rather sensitive for differences of small areas. Besides a shift 
of A,B, to the right or to the left enlarges one area and diminishes 
the other so that even a slight deviation from the correct position 
makes itself felt. In the same way the next step B,CiC2 is 
drawn with its vertical part B,C; in such a position that the 
areas on both sides are equal. The integral of the given curve 
as far as K, will again have the same value as that of the stepping 
line as far as K,. And so on for the other steps. The integral 
of the stepping line is constructed in the way shown. It is 
represented by a broken line beginning at the foot of the ordinate 
of A,;. The corners lie on the vertical parts of the steps or 
their prolongations. It is readily seen that the broken line con- 
sists of a series of tangents of the integral curve 


y= [ ferdr! 


and that their points of contact with the integral curve lie on 
the same verticals as the points A1, Kz, K., ete. (In Fig. 73 these 
points are denoted 0, 2, 3, ---.) That these points lie on the 
integral curve follows from the arrangement of the steps which 
make the integral of the given function at K;, K., --+ equal tothe 
integral of the stepping line. Now in the points Aj, Kz, K, --- 
the ordinates of the given curve coincide with those of the 
stepping line. Hence both integral lines must for these abscissas 
have the same direction. 


1In Fig. 73 the lower limit is 0. 
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Having constructed the broken line and marked the points 
2, 3, 4, --+ (Fig. 73), the integral curve is drawn with a curved 
; ruler so as to touch the broken line in the polntsy 0,22, 85 se 
{ As the given curve does not change its ordinate abruptly the 
integral curve does not change its direction abruptly. The 
drawing shows how well the integral curve is determined by the 
broken line. There is practically no choice in drawing it any 
other way without violating the conditions. 
The ordinate of the integral curve is measured in the same 
: unit as the ordinate of the given curve y = f(z). It may some- 
times be convenient to draw the ordinates of the integral curve 
in a scale different from that of the ordinates of the given curve. 
For instance the value of the integral may become so large that 
measured in the same unit the ordinates of the integral curve 
would pass the boundaries of the drawing board, or else they may 
be so small that their changes cannot be measured with sufficient 
accuracy. In the first case the scale is diminished, in the latter 
case it is enlarged. This is done by altering the position of the 
point P, the center of the pencil of rays that define the directions 
of the broken line. If P approaches O the directions Pa, Pf, --- 
become steeper to the same degree as if keeping P unchanged we 
had increased the ordinates of A142, B, Bo, - ++ in the inverse pro- 
- portion of the two distances PO. Hence by diminishing the 
- distance PO the ordinates of the resulting broken line are enlarged 
in the inverse proportion. On the other hand, by increasing the 
distance PO the ordinates of the resulting broken line are di- 
minished in the inverse proportion of the distances, because the 
change of the directions Pa, PB, --- caused by a longer distance 
PO is the same as if the ordinates of A,Ao, B,B:, --- were di-' 
minished in the inverse proportion. The broken line constructed 
by means of the longer distance P’O will therefore be the same as 
if the ordinates of the stepping line were diminished. It therefore 
leads to an integral curve whose ordinates are diminished in the 
same proportion (Fig. 74). 
The graphical integration of 
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Y= [sede 


is not limited to values z > a. The method is just as well applic- 
able to the continuation of the integral curve for <a. The 


Fia. 74. 


steps have only to be drawn from right to left. The lower limit 
a determines the point where the integral curve intersects the 
axis of x. 

There is a method for the construction of the vertical parts 
of the steps, which may in some cases be useful, though as a rule 


we may dispense with it and fix their position by estimation. 


Suppose that A and B (Fig. 75) 


are two points where the curve is 
intersected by the horizontal parts 
of two consecutive steps and that 
the curve between A and B is a 
parabola whose axis is parallel to 
the axis of x. The position of the 
vertical part of the step between 4 
and B can be then found bya simple 
construction. Through the center C of the chord AB (Fig. 75) 
draw a parallel CD to the axis of x, D being the point of inter- 
section with the parabola. The vertical part EH of the step in- 
tersects CD ina point whose distance from C is twice the distance 


| 
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from D. That this is the right position of EH is shown as soon 
as we can prove that the area ADBGA is equal to the rectangle 
EHBG. The area ADGBA can be divided in two parts, the tri- 
angle ABG and the part ADBCA between the curve and the 
chord. The triangle is equal to the rectangle FIBG, while ADBCA 
is equal to two thirds of the parallelogram MNBA, and hence 
equal to the rectangle HHIF. Both together are therefore equal 
to the rectangle EH BG, and the two areas between the stepping 
line and the curve on both sides of EH are thus equal. 

If the curve between A and B is sup- 
posed to be a parabola with its axis par- 
allel to the axis of ordinates the con- 
struction has to be modified a little. 
Through the center C of the chord AB 
(Fig. 76) draw a vertical line CD as far 
as the parabola. On CD find the point 
K whose distance from C is double the 
distance from D and draw through it a 
parallel to the chord AB. This parallel 
intersects a horizontal line through C at a point LZ. Then EH 
must pass through L. This may be shown in the following way. 
The area between the parabola ADB and the chord AB is equal 
to two thirds of the parallelogram MNBA, MN being the tan- 
gent to the parabola at the point D. If D’ is the point of inter- 
section of VN and the horizontal line through C, we have evi- 
dently 


CL = 2CD’. 


Therefore the rectangle EHIF is equal to the area ADBA be- 
tween the parabola and the chord and EHBG is equal to ADGBA. 

Any part of a curve can be approximated by the arc of a 
parabola with sufficient accuracy if the part to be approximated 
is sufficiently small. When the direction of the curve is nowhere 
parallel to the axis of codrdinates, both kinds of parabolas may 
be used for approximation, those whose axes are parallel to the 
axis of a and those whose axes are parallel to the axisof y. But 


i 
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when the direction in one of the points is horizontal (Fig. 76), 
we can only use those with vertical axes and when the direction 
in one of the points is vertical we can only use those with hori- 
zontal axes. Accordingly we have to use either of the two con- 


structions to find the position of the vertical part of the step. — 


Do not draw your steps too small. For, although the difference 
between the broken line and the integral curve becomes smaller, 
the drawing is liable to an accumulation of small errors owing 
to the considerable number 
of corners of the broken 
line and little errors of 
drawing committed at the 
corners. Only practical ex- 
perience enables one to find 
the size best adapted to 
the method. 

Fic. 77. Statical moments of areas 

may be found by a double 

graphical integration. Let us consider the area between the curve 

y = f(x) (Fig. 77), the axis of x and the ordinates corresponding 

tow = Oand «= £ The statical moment with respect to the 

vertical through 2 = & is the integral of the products of each 
element ydz and its distance § — x from the vertical 


é 
M= f (E — x)yda. 


Let us regard M as a function of £ and differentiate it: 


dM od 
= (E— eet f eG ayde 


= 0+ [ ydx. 


That is to say, a graphical integration of the curve y = f(a) 
beginning at 2 = 0 furnishes the curve whose ordinate is 

aM 

dé * 


r 
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Hence a second integration of this latter curve will furnish 
the curve M as a function of ¢& As M vanishes for & = 0 the 
second integration must also begin at the abscissa 2 = 0. 


Fig. 78 shows an example. Each ordinate of the curve found 
by the second integration is the statical moment of the area on 
the left side of it with respect to the vertical through this same 
ordinate. The ordinate furthest to the right is the statical 
moment of the whole area with respect to the vertical on the 
right. The statical moment of the whole area with respect to a 
vertical line through any point 2; is the integral 


f ; (a, — x)ydz. 


Considered as a function of 2, its differential coefficient is 


Ey f 
f pes (a, — x)ydx 4, ydx. 


That is to say, the differential coefficient is independent of 1, 
hence the statical moment is represented by a straight line. As 
its differential coefficient is represented by a horizontal line 
through the last point on the right of the curve 


& 
[ ydz, 


‘ 
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the direction of the straight line is found by drawing a line 
through P and through the point of intersection Q of the hori- 
zontal line and the axis of ordinates (Fig. 78). The position of 
the straight line is then determined by the condition that 


{ om — x)ydx 


for x; = £ is equal to the statical moment 
r€ 
M@® = (— ayde, 


We have therefore only to draw a parallel to PQ through the 
last point R of the curve for M(é) found by the second integration. 
The ordinates of this straight line for any abscissa 2; represent 
the values of 


7 
: i] (a1 — x)ydx 


measured in the unit of length of the ordinates. The point of 
intersection E with the axis of x determines the position of the 
vertical in regard to which the statical moment is zero, that is to 
say, the vertical through the center of gravity. 

The moment of inertia of the area 


& 
{i ydx 
0 


about the axis « = € is found ina similar way. It is expressed 
by the integral 


F 
T= i (§ — wx)*ydz. 


Considered as a function of £ we find by differentiation 


dT = d 
aa [Ce Pyle t+ f di (£ — x)2ydx 


é 
= 0+2[ aude. 


os 
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That is to say, the differential coefficient is equal to double the 


statical moment about the same axis. This holds for every value 
of £ Hence we obtain $7 as a function of £ by integrating 
the curve for M(£). For & = 0 we have JT = 0, so that the curve 
begins on the axis of 2 at & = 0. ; 


The integral 
J ydx 


is zero for «=a. The curve representing the integral has to 
intersect the axis of x at x = a (admitting values of z > a and 
x <a), and it is there that we begin the construction of the 
broken line. If instead we begin it at the point 2 = a, y = ¢, 
the only difference is that the whole integral curve is shifted 
parallel to the axis of ordinates by an amount equal to c upwards 
if c is positive, downwards if it is negative. But the form of the 
curve remains the same. It is different when this curve is 
integrated a second time. For instead of 


fi ydx 
ii ydx + ¢. 


The ordinate of the integral curve is therefore changed by an 
amount equal to c(a — a) and besides if the second integral curve 
is begun at x = a, y= G instead of =a, y= 0 the change 
amounts to 


we now integrate 


c(x a a) + C1, 
so that the difference between the ordinates of the new integral 
curve and the ordinates of the straight line 
y=ex—ayta 


is equal to the ordinates of the first integral curve (Fig. 79). 
This effect of adding a linear function to the ordinates of the 
integral curve is also attained by shifting the pole P upward or 
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downward. For it evidently comes to the same thing whether 
the curve to be integrated is shifted upward by the amount ¢ or 
whether the point P is moved downward by the same amount, so 
that the relative position of P and the curve to be integrated 
is the same as befote. Changing the ordinate of P by — cadds — 


cn 
(a) dada 
ia 


KK 
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c(x — a) to the ordinates of the integral curve. c(2 — a) is the 
ordinate of a straight line parallel to the straight line from the 
new position of P to the origin. 

By this device of shifting the position of P upward or down- ~ 
ward the integral curve may sometimes be kept within the 
boundaries of the drawing without any reduction of the scale of 
ordinates. A good rule is to choose the ordinate of P about 
equal to the mean ordinate of the curve to be integrated. The 
ordinates of the integral curve will then be nearly the same at 
both ends. The value of the integral 


Hi ydx 


is equal to the difference between the ordinates of the integral 
curve and the ordinates of a straight line parallel to PO through 
the point of the integral curve whose abscissa is a. 
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When the ordinate of P is accurately equal to the mean ordinate 
of the curve to be integrated for the interval x = a to b the 
ordinates of the integral curve will be accurately the same at the 
two ends. But we do not know the mean ordinate before having 
integrated the curve. 

After having integrated we find the mean ordinate for the 
interval 2 = a to b by drawing a straight line through P parallel 
' to the chord AB of the integral curve, A and B belonging to the 
abscissas e=a and z=). This 
_line intersects the axis of ordi- 
nates at a point whose ordinate 


is the mean ordinate. A 4 iB 
Suppose a beam AB is sup- i ‘ [ 
ported at both ends and loaded Fic. 80. 


by a load distributed over the 
beam as indicated by Fig. 80. That is to say, the load on dz is 
measured by the area ydz. Let us integrate this curve graph- 
ically, beginning at the point A with P on the line AB. The 
final ordinate at B 

f ydx 


gives the whole load and is therefore equal to the sum of the two 
reactions at A and B that equilibrate the load.. Integrating this 
curve again we obtain the curve whose ordinate is equal to 


E 
i Vas: 
if ydz 


The ordinate of this curve at any point x = & represents the 

statical moment of the load between the verticals x = a and 

x = £ about the axis x = &. Its final ordinate BM, Fig. 81,.1s’ 

the moment of the whole load about the point B, and as the reac- 

tions equilibrate the load it must be equal to the moment of the 
9 


Y being written for 
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reactions about the same point and therefore opposite to the 


moment of the reaction at A about B. If the reaction at A is 
denoted by F, we therefore have 


b 
Fab — a) = [Yate 


That is to say, F, is equal to the mean ordinate of the curve 


é 
Y= | ydz 


in the interval = a to b. The mean ordinate is found by 
drawing a parallel to AM through P which intersects the vertical 
through A at the point F so that AF = F,. As DB is equal to 


| ds i 
SSF (ee) decette 
| a a; : M 


= 
~ 
oR 
~. 
ee 


1 
! 1 
ea 
eee ee ee 
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the sum of the two reactions a horizontal line through F will 
divide BD into the two parts BG = F, and GD = F,. 


Shifting the position of P to P’ on the horizontal line FG 
and repeating the integration é 


; 
f Ydxz, 


we obtain a curve with equal ordinates at both ends. If we 
begin at 4 it must end in B. Its ordinates are equal to the 
difference between the ordinates of the chord 4M and the curve 
AM (Fig. 81), and represent the moment about any point of 


ry Tre 0 i 
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the beam of all the forces on one side of the point (load and 
reaction). 

The area of a closed curve may be found by integrating over the 
whole boundary. Suppose z = a and x = b to be the limits of 
the abscissas of the closed curve, the vertical x = a touching the 
curve at A and the vertical z = 6 at B (Fig. 82). By A and B 
the closed curve is cut in two, both parts connecting A and B. 
Let us denote the upper part by y = f,(x) and the lower part 
by y = f2(x). The whole area is then equal to the difference 


x 
JFa(@dee 


Fie. 82. 


[n@d—[ pepae, 


or equal to , : 
i) fila)dx + if fo(x)dx. 


We begin the integral curve over the upper part at the vertical 
x2 = aata point £, the ordinate of which isarbitrary, and draw 
the broken line as far as F on the vertical « = b (Fig. 82). Then 
we integrate back again over the lower part, continuing the 
broken line from F to G. The line HG measured in the unit of 
length set down for the ordinates is equal to the area measured 
in units of area, this unit being a rectangle formed by PO and 
the unit of ordinates. That is to say, the area is equal to the 
area of a rectangle whose sides are PO and EG. 
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The method is not limited to the case drawn in Fig. 82, where 
the closed curve intersects any vertical not more than twice. A 
more complicated case is shown in Fig. 83. But in all those cases 


Fic. 83. 


where the object is not to find the integral curve but only to find 
the value of the last ordinate the method, cannot claim to be 
of much use, because it cannot compete with the planimeter. 


For the construction of the broken line we have drawn the 
steps in such a manner that the areas on both sides of the vertical 
part of a step between the curve and the stepping line are equal. 
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It would have also been admissible to construct the stepping 
line in such a way that the areas on both sides of the horizontal 
part of a step are equal (Fig. 84). Only the broken line would 
consist of a series of chords instead of a series of tangents of the 
integral curve. The points K., Ks, ---, where the horizontal 
parts of the steps intersect the curve would determine the ab- 
scissas of the points of the integral curve, where its direction is 
parallel to the direction of the broken line. But this forms very 
little help for drawing the integral curve. That is the reason 
why the former method where the broken line consists of a series 
of tangents is to be preferred. However where the object is only 
to find the last ordinate of the integral curve the two methods 
are equivalent. 

§ 14. Graphical Differentiation.—The graphical differentiation 
of a function represented by a curve is not so satisfactory as the 
graphical integration because 
the values of the differential 
coefficient are generally not 
very well defined by the curve. 
The operation consists in 
drawing tangents to the given 
curve and drawing parallels 
through P to the tangents 
(Fig. 85). The points of in- 
tersection of these parallels Fra. 85. 
with the axis of ordinates fur- 
nish the ordinates of the curve representing the derivative. 
The abscissa to each ordinate coincides with the abscissa of the 
point of contact of the corresponding tangent. The principal 
difficulty is to draw the tangent correctly. As a rule it can be 
recommended to draw a tangent of a given direction and then 
mark its point of contact instead of trying to draw the tangent 
for a given point of contact. A method of finding the point of 
contact more accurately than by mere inspection consists in 
drawing a number of chords parallel to the tangent and to 
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bisect them. The points of bisection form a curve that inter- 


sects the given curve at the point of contact (Fig. 86). When a _ 


number of tangents are drawn, their points of contact marked 
and the points representing the differential coefficient constructed, 
the derivative curve has to be 
drawn through these points. 
This may be done more accur- 
ately by means of the stepping 
line. The horizontal parts of 
the steps pass through the 
points while the vertical parts 
lie in the same vertical as the 
point of intersection of two 
consecutive tangents. The derivative curve connects the points 
in such a way that the areas between it and the stepping line are 
equal on both sides of the vertical parts of each step. Thus 
the result of the graphical differentiation is exactly the same 


Fic. 86. 


Fic. 87. 


figure that we get by integration, only the operations are carried 
out in the inverse order. 

A change of the distance PO (Fig. 87) changes the ordinates 
of the derivative curve in the same proportion and for the same 
reason that it changes the ordinates of the integral curve when we 


uth Ahan, nner lee teal 


— 
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are integrating, but in the inverse ratio. Any change of the or- 
dinate of P only shifts the curve up or down by an equal amount, 
so that if we at the same time change the axis of x and draw it 
through the new position of P the ordinates of the curve will 
remain the same and will represent the differential coefficient. 

When a function f(x, y) of two variables is given by a diagram 
showing the curves f(z, y) = const. for equidistant values of 
f(x, y) the partial differential coefficients can be found at any 
point 2, yo by means of drawing curves whose ordinates represent 
f(x, yo) to the abscissa x or f(xo, y) to the abscissa y and applying 
the methods explained above. For this purpose a parallel is 
drawn to the axis of 2, for instance, through the point 2, yo 
and at the points where it intersects the curves f(z, y) = const. 
ordinates are erected representing the values of f(z, yo) in any 
convenient scale. A smooth curve is then drawn though the 
points so found and the tangent of the curve at the point 2» 
furnishes the differential coefficient df/dx for x = a, y = Yo. 

The differential coefficients df/dx, af/dy are best represented 
graphically by a straight line starting from the point 2, y to 
which the differential coefficients correspond, and of such length 
and direction that its orthogonal projections on the axis of x 
and y are equal to df/dzx and df/dy. This line represents the 
gradient of the function f(z, y) at the point a, y. It is normal 
to the curve f(z, y) = const. that passes through the point 2, y, 
its direction being the direction of steepest ascent. Its length 
measures the slope of the surface z = f(z, y) in the direction of 
steepest ascent. This is shown by considering the slope in any 
other direction. Let us change x and y by 


rcosa, rsina 
and consider the corresponding change 
Az = f(«+rcosa,y+rsina).— f(x, y) 
of he function. By Taylor’s theorem we can write it 


1 See Chap. IT, § 10. 
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af of 
og f cos a+ ay 


a is the direction from the point z, y to the new point z + r cos a, 
y+ rsina and r is the distance of the two points. Dividing 
Az by r and letting r approach to zero we find 


r sin a + terms of higher order in r, 


This expression measures the slope of the surface z = f(xy) 
in the direction a. Now let us introduce the length / and the 
angle \ of the gradient, and write 

of 


of ; 
es l cos X, af I sin X. 


Then we have 


fe) Ofivaxs 
oF cos a + 7 sine = 1 cos (a — X). 
That is to say, the slope in any direction a@ is proportional to 
cos (a — \), it is a maximum in the direction of the gradient 
(a = \) and zero in a direction perpendicular to it and negative 
in all directions that form an obtuse angle with it. When all 


three codrdinates are measured in the same unit, the length of 


J measured in this unit is equal to the tangent of the angle of — 


steepest ascent. Hence the length of the gradient varies with 
the unit of length. When ‘the unit of length in which the values 
of f(xy) are plotted is kept unaltered, while we change the unit 
of length corresponding to the values 2 and y, the length of the 
gradient varies with the square of the unit of length. 

$15. Differential Equations of the First Order.—In the problem 
of solving a differential equation of the first order 


d 
72 @y) 


by graphical methods the first question is how to represent 
the differential equation graphically. If 2 and y are meant to 
be the values of rectangular codrdinates, the geometrical meaning 


Ee 


a ee 
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of the differential equation is that at every point x, y, where 
f(x, y) is defined, the equation prescribes a certain direction for — 
the curve that satisfies it. Let us suppose curves drawn through 
all those points for which f(a, y) has certain constant values. 
Each curve then corresponds to a certain direction or the opposite 
direction. Let us distinguish the curves by different numbers or 
letters and let us draw a pencil of rays together with the curves 
and mark the rays with the same numbers or letters in such a way 
that each of them shows the direction corresponding to the 


2nd approximation 
(Integration of the 
curve below) 
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curve marked with that particular number or letter (Fig. 88). 
Our drawing of course only comprises a certain region in which 
we propose to find the curves satisfying the differential equation. 
It may be that f(zy) is defined beyond the boundaries of our 
drawing. Those regions have to be dealt with separately. 

The graphical representation of the differential equation in 
the region considered consists in the correspondence between 
the curves and the rays. It is important to observe that this 
representation is independent of the system of coérdinates by 
means of which we have deduced the curves from the equation 


oY = Fay). 


: 
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We can now introduce any system of coérdinates £ 7 and find : 
from our drawing the equation 


d 
- = ¢(én), 


that is to say, we can find the value of ¢(£, 7) at any point & 7 
of our drawing. If, for instance, the unit of length is the same 
for £ and 7 we draw a line through the center of the pencil of rays 
in the direction of the positive axis of £ and a line perpendicular to 
it at the distance 1 from the center. The segment on the second 
line between the first line and the point of intersection with one 
of the rays measured in units of length and counted positive in 
the direction of positive 7 furnishes the value of ¢(&, 7) for all 
the points £, » corresponding to that particular ray. In this 
respect the graphical representation of a differential equation 
is superior to the analytical form, in which certain coérdinates 
are used and the transformation to another system of codrdinates 
requires a certain amount of calculation. 

Now let us try to find the curve through a given point P on 
the curve marked (a) (Fig. 88) that satisfies the differential equa- 
tion. We begin by drawing a series of tangents of a curve 
that is meant to be a first approximation. Through P we draw 
a parallel to the ray (a) as far as the point Q somewhere in the 
middle between the curves (a) and (b). Through Q we draw a 
parallel to the ray (b) as far as R somewhere in the middle 
between the curves (b) and (c). Through R we again draw a 
parallel to the ray (c) and so on. The curve touching this 
broken line at the points of intersection with the curves (a), 
(b), -:+ is a first approximation. But we need not draw this 
curve. In order to find a better approximation we introduce a 
rectangular system of codrdinates 2, y, laying the axis of 2 some- 
what in the mean direction of the broken line. Let us denote 
by ym the function of x that corresponds to the curve forming the 
first approximation. The second approximation y is then ob- 
tained as an integral curve of f(x, y:), that is, of dy,/dx 


fal bes a 
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n= ot fie, war, 


denoting by 2p, yp, the codrdinates of P. For this purpose the 
curve whose ordinates are equal to f(z, y:) or dy;/dx has to be con- 
structed first. The values of f(z, y:) are found immediately at 
the points where the first approximation intersects the curve 
(a), (b) --- by differentiation in the way described above. A 
line is drawn through the center of the pencil of rays parallel to 
the axis of x and a line perpendicular to it at a convenient dis- 
tance from the center. This distance is chosen as the unit of 
length. The points of intersection of this line with the rays de- 
termine segments whose lengths are equal to the values of f(a, y1) 
on the corresponding curves. These values are plotted as ordi- 
nates to the abscissas of the points where the first approximation 
intersects the curves (a), (b), --- and a curve 


Y= J; YM) 


is drawn (Fig. 88). This curve is integrated graphically begin- 
ning at the point P and the integral curve is a second approxi- 
mation. Again we need not draw the curve. The broken line 
suffices, if we intend to construct a third approximation. In 
this case we have to repeat the foregoing operation. ‘This can 
now be performed much quicker than in the first case because the 


values of f(z, y) on the curves (a), (b), --+ have already been 
constructed and are at our disposal. In order to find the curve 
Y= i (x, ye) 


we have only to shift the same ordinates to new abscissas and 
make these coincide with the abscissas of the points where the 
second approximation intersects the curves (a), (b), ---. The 


curve 
Y= jae cP) 


is then drawn and integrated graphically, beginning at the point 
leg 
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Suppose now the integral curve did not differ from the second 
approximation, it would mean that 


nat { fewde 
or that 
d 
ae = 1 w), 


that is to say, that y satisfies the differential equation. 

If there is a perceptible difference the integral curve represents 
a third approximation. It has been shown by Picard that pro- 
ceeding in this way we find the approximations (under a certain 
condition to be discussed presently) converging to the true solu- 
tion of the differential equation, so that after a certain number 
of operations the error of the approximation must become 
imperceptible. 

Denoting by y, the function of the nth approximation we have 


Yr = Yo =e {i i E> Bi Nee: 


The true solution with the same initial conditions y = yy for 
x = x, satisfies the equation 


Y = Yp + | F(a, y)dz. 


Hence 


yoy = [Li yn) — fle, Wilde, 


* F(x, Yn) — F&, 
= Yn — Y 


Let us now suppose that the absolute value of 


Yn— ¥ 


or 


’ 


for all the values of a, y, yn within the considered region does 


—" 


le 
‘n 
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not surpass a certain limit M, then it follows that a certain relation 
- must exist between the maximum error of y,, which we denote by 
en and the maximum error of yn41, which we denote by en41. 
The absolute value of the integral not being larger than 


Me, |x — 2, | 
(|x — 2, | denoting the absolute value of x — z,) we have 
Catt — | 2 a, | ey, 


Hence as long as the distance 2 — 2, over which the integration 
is performed is so small that 


M\|x—2|<k <1, 


k being a constant smaller than one, the error of yn41 cannot be 
larger than a certain fraction of the maximum error of yp. 
But in the same way it follows that the error of y, cannot be 
larger than the same fraction of the maximum error of yp_1, and 


so on, so that 
Oo os Neg = he, 1 > = Ee, 


But as e is a constant and fk a constant smaller than one, k”e 

must be as small as we please for a sufficient large value of n. 

That is to say, the approximations converge to the true solution. 
M being a given constant the condition of convergence 


M\|x—2,|<k<1 


limits the extent of our integration in the direction of the axis of z. 
But it does not limit our progress. From any point P’ that we 
have reached with sufficient accuracy we can make a fresh start, 
choosing a new axis of x suited to the new situation. As a 
rule it does not pay to trouble about the value of M and to try 
to find the extent of the convergence by the help of this value. 
The actual construction of the approximations will show clearly 
enough how far to extend the integration. As far as two consecu- 
tive approximations show no difference they represent the true 


curve. 


7 
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Suppose that ; 
f(a, Yn) — f(a, y) ; 
Yn—Y 


has the same sign for all values x, y, yn concerned. Say it is 
negative. Suppose further that y, — y is of the same sign for 
the whole extent of the integration 


Ynti— Y= i ie = — i 2 (yn — y)da; 


that is to say, the approximative curve Y» is all on one side of the 
true curve. Then if x — 2p is positive, ¥ni1 — y must evidently 
be of the opposite sign from y, — y, or the approximative curve 
Yn+1 1s all on the other side of the true curve from y,. For these 
and all following approximations the true curve must lie between 
two consecutive approximations. If the first approximation y; is 
all on one side of the true curve the theorem holds for any two 
consecutive approximations. This is very convenient for the esti- 
mation of the error. 
In Fig. 88 
F(X, Yn) — F(, y) 
Yn— Y 


is negative from the point P as far as somewhere near S. The 
first approximation is all on the upper side of the true curve. 
Therefore the second approximation must be below the true 
curve at least as far as somewhere near S. 

When the sign is positive the same theorem holds for negative 
values of e — x». If the integration has been performed in the 
positive direction of x, it may be a good plan to check the result 
by integrating backwards, starting from a point that has been 
reached and to try if the curve gets back to the first starting 
point. In this direction we profit from the advantage of the 
true curve lying between consecutive approximations and are 
better able to estimate the accuracy of our drawing. 

We have seen that the convergence depends on the maximum 


aw 
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absolute value of 
f(z, Yn) ze f(2, y) 
Yn— Y 

for all values of x, y, yn concerned. In order to find the maximum 
value we may as well consider 

of 

ay 
for all values of x, y within the region considered. For if we 
assume df/dy to be a continuous function of y, it follows that 
the quotient of differences 


F(X, Yn) — F(@, Y) 
Yrn— ¥ 

must be equal to of/dy taken for the same value of z and a value 
of y between y and y,. This is immediately seen by plotting 
f(x, y) as ordinate to the abscissa y for a fixed value of xz. The 
value of the quotient of differences is determined by the slope 
of the chord between the two points of abscissas y and yn. The 
slope of the chord is equal to the slope of the curve at a certain 
point between the ends of the chord. The value of df/dy at this 
point is equal to the value of 


F(&, yn) — F(a, Y) 
i ae 

Now let us consider how the coérdinate system may be chosen 
in order to make df/dy as small as possible and thus obtain the 
best convergence. For this purpose let us investigate how the 
value of df/dy changes at a certain point, when the system of 
codrdinates is changed. 

Let us start with a given system of rectangular coérdinates &, 
n with which the differential equation is written 


d 
7: = o(&, 7). 


The direction of the curve satisfying the differential equation 
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: 


4 
s 
; 


forms a certain angle a with the positive axis of & determined by 


dn 
1S oe o(, n) 
(assuming the codrdinates to be measured in the same unit). 
Now let us introduce a new system of rectangular codrdinates — 
f 


x, y connected with the system £, 7 by the equations | 
x= ~Ecoswt+ sina, 
y = —ésinw+ 7 Cosa, 
which are equivalent to 
&= rcosw — ysina, 
n= xsinw+ y cosa, 


w being the angle between the positive direction of x and the 
positive direction of £, counted from & towards 2 in the usual way. 

The angle formed by the direction of the curve with the positive 
direction of the axis of x is a — w, and therefore 


oY = te (a—«) =f »). 


Consequently we obtain for a given value of w 


of 1 0a 


dy cos? (a — w) dy’ 


or remembering that a is given as a function of £ and 7, 


of 1 da , da 
ay = waza (-% sin w + se cose ), 
For simplicity’s sake we shall assume that the axis of € is the 
tangent of the curve ¢(é, 7) = const. that passes through the 
given point, so that da/dé = 0. 
We then have 
of re 1 Oa 


dy cos? (a — w) dn pee 


and our object is to find how éf/dy varies for different values of 


in the direction of the gradient. 


; outs oe co it us nes ie ; 


Dp wae 


a, which we represent by a straight line drawn 


rigin A (Fig. 89) perpendicular to the curve a = const. 
= const. 


wae 
0 restriction to assume the value of da/dn positive; it 


ns that the direction of the positive axis of 7 is chosen 


\ 
\ 
\ 
\ 
\ 


\ 
\ 
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Let us draw the line AB (Fig. 


89) in the direction of the peotys axis of £ and of the same length 
_as the gradient. 


In order to show the values of of/dy for the different positions 
of the axis of 2 let us lay off the value of daf/dy as an abscissa. 
For instance for w = a, of/dy assumes the value 


ae os 
eal a. 
On 


- The abscissa corresponding to this value is AB’ (Fig. 89), the 


10 
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orthogonal projection of AB on the axis of . For any other 
position AC (Fig. 89) corresponding to some other value of o, 
we find da/dn cos w by orthogonal projection of AB on AC. Then 
the division by cos (a — w) furnishes AC’ and a second division 
by cos (a — w) leads to AC. Thus a certain curve can be 
constructed whose polar codrdinates are r = Of/dy and w, the 
equation in polar coérdinates being : 


——————eEoo 


da COS @ 
On cos? (a — a) 


to) 
o— or [rcos (a — w)P= = COS w. 


In rectangular coérdinates £, 7 the equation assumes the form 
0a 


(cos aé + sin an)? = rs g, 


This shows that the equation is a parabola, the axis of which is 
perpendicular to the direction a. AB’ isa chord and the gradient 


4 
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A(risa tangent of the parabola. Bisecting AB’ in E, drawing EK 
perpendicular to AB’ as far as the axis of n and bisecting EK in 
D, we find D the apex of the parabola. The three points A, B’, D 
together with the gradient will suffice to give us an idea of the 


size and sign of df/dy for the different positions of the positive 
axis of 2. 
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df/dy vanishes when the axis of a is perpendicular to the curve 
a = const., so that it seems as if this were the most favorable 
position. We must, however, bear in mind that the axis‘of a 
is kept unaltered for a certain interval of integration. When we 
pass on to other points the axis of x is no longer perpendicular 
to the curve a = const. there. The position of the axis of z is 
good when the average value of df/dy is small. In Fig. 90 the 
parabolas are constructed for a number of points on the first 
approximation of a curve satisfying the differential equation. 

If we want to make use of the parabolas to give us the numerical 
values of df/dy the unit of length must also be marked in which 
the codrdinates are measured. The numerical value of df/dy 
varies as the unit of length and therefore the length of the line 
representing it must vary as the square of the unit of length. 
But if we draw a line whose length measured in the same unit is 


r ae this line would be independent of the unit of 


length. For if ] is the line representing the unit of length and 


equal to 


: 1 
I’, I’ the lines representing the values af/dy and aflay’ of/ay 


would be the ratio /’/] and =~ the ratio l’’/l; hence l’” = P/I’. 


a 
Since J’ varies as 2 with the change of the unit of length /’’ is 
independent of the unit of length. This line /’’ represents the 
limit beyond which the product 


of jy’ 

oy 
becomes greater than 1. If af/dy remained the same this would 
mean the limit beyond which the convergence of the process of 


; 1 
approximation ceases. We might lay off the length of afley in 


the different directions in the same way as df/dy has been laid 
off. The result is a curve corresponding, point by point, to the 
parabola, the image of the parabola according to the relation of 
reciprocal radii. But all these preparations as a rule would not 
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pay. It is better to attack the integration at once with an axis _ 


of « somewhat perpendicular to the curves a = const. as long 
as the direction of the curve forms a considerable angle with 
the curve a = const. and to lose no time in troubling about the 
very best position. The convergence will show itself, when the 
operations are carried out. When the angle between the direction 
of the curve that satisfies the differential equation and the curve 
a = const. becomes small the apex of the parabola moves far 
away and when the direction coincides with that of the curve 
a = const. the parabola degenerates into two parallel lines per- 
pendicular to the direction of the curve a = const. In this case 
the best position for the axis of z is in the direction of the curve 
a = const. Without going into any detailed investigation about 
the best position of the axis of x we can establish the general rule 
not to make the axis of x perpendicular to the direction of the 
curve satisfying the differential equation, that is to say, not to 
make it parallel to the axis of the parabola. But we hardly need 
pronounce this rule. In practice it would enforce its own observ- 
ance, because for that position of the axis of 2 not only of/dy but 
also f(z, y) are infinite and it would become impossible to plot 
the curve Y = f(a, 4). 


There is another graphical method of integrating a differential _ 


equation of the first order 


which in some cases may well compete with the first method. 
Like the first it is the analogue of a certain numerical method. 

The numerical method starts from given values 2, y and cal- 
culates the change of y corresponding to a certain small change 
of z. Let h be the change of x and k the change of y, so that 
a+ h,y + kare the coérdinates of a point on the curve satisfying 
the differential equation and passing through the point 2, y. kas 
calculated in the following manner. We calculate in succession 
four values hy, ko, ks, ks by the following equations— 


. 
7 
é 
; 
. 
i 
4 


Ae ‘— 
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ky a A (a, y) h, 
Brn, ts 
ky = i(e+$.y4 aE 


= h hy 
kg = i(2 + yt ae 
kg = f(a th, y + ks)h. 
We then form the arithmetical means 


hk. + ks ky + kg 
Ne as and gate a 


and find with a high degree of approximation as long as h is 
not too large 
k= p+¥a-2)s 
The new values 
A=e+h, Y=y+tk 

are then substituted for x and y and in the same way the codrdi- 
nates of a third point are calculated and so on. 

This calculation may be performed graphically in a profitable 
manner, if the function f(z, y) is represented in a way suited to 


ym 
Fig. 91. 
1See W. Kutta, Zeitschrift fiir Mathematik und Physik, Vol. 46, p. 443. 


the purpose. Let us suppose a number of equidistant parallels 
to the axis of ordinates: x = %, ® = %1, © = %, X= XH, °** 
Along these lines f(z, y) is a function of y. Let us lay off the — 
values of f(x, y) as ordinates to the abscissa y, the axis of y being ‘ 
taken as the axis of abscissas. We thus obtain a number of 
curves representing the functions f(2o, y), f(a, y), f(a, Y), ++" 
Starting from a point A(zo, yo) on the first vertical 2 = ao (Fig. 
91) we proceed to a point B; on the vertical 2 = 2 in the following 
way. By drawing a horizontal line through A we find the 
point A’ on the curve representing f(ao, y). Its ordinate is equal 
to f(xo, yo). Projecting the point A’ onto the axis of x we find A” 
and draw the line PA”. P is a point on the negative side of — 
the y-axis and PO is equal to the unit of length by which the 
lines representing f(x, y) are measured. Thus 


OA"/PO = f(xo, Yo). 
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Now we draw AB, perpendicular to PA”, so that if h and ky 
denote the differences of the codrdinates of A and B, we have 


ky/h = OA" /PO, 


ky = f (xo, Yoh. 


From C; the point of intersection of the line 4B, and the vertical 
x = x, we find C,/ and C,;” in the same way as we found A’ and A” 
from A, only that Cj’ is taken in the curve representing the 
values of f(a1, y), and draw the line AB, perpendicular to PO”. 
Denoting the difference of the ordinates of A and B; by ky we have 


Is “OC,” mee 


or 
h k 
ky = ic + 5? Yo + 3)h 


From C, the point of intersection of the line AB, and the 
vertical 2 = 2; we find in the same way a point B; on the vertical 
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x = a and the difference k3; between the ordinate of B; and that 
of A is 


h k. 
ks =s(2 +5) Yo ate oh 


From B3 we pass horizontally to B;’ on the curve representing 
f(ay) and vertically down to B;’’.. The line AB, is then drawn 
perpendicular to PB;’’, so that the difference ky between the 
ordinates of B, and A is 


ky = f(to + hy y + hes)h. 


The bisection of B,B; and of B,B, gives us the points FE; and Ly 
and the point B is taken between FE; and Ey, so that its distance 
from F) is half its distance from E,. The point B is with a high 
degree of approximation a point of the curve that passes through 
A and satisfies the differential equation. 

B is then taken as a new point of departure instead of A, and 
in this manner a series of points of the curve are found. 

In order to get an idea of the accuracy attained the distance 
of the vertical lines is altered. For instance, we may leave out 
the verticals x = 2, and x = 23, and reach the point on the 
vertical x = x, in one step instead of two. The error of this 
point should then be about sixteen times as large as the error on 
the same vertical reached by two steps, so that the error of the 
latter should be about one-fifteenth of the distance of the two. 
If their distance is not appreciable the smaller steps are evidently 
unnecessarily small. 

The values of f(z, y) may become so large that an incon- 
veniently small unit of length must be applied to plot them. In 
this case 2 and y have to change parts and the differential 
equation is written in the form 


dx 1 
dy fy)” 


The values of 1/f(z, y) are then plotted for equidistant values of 
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y as ordinates to the abscissa z and the constructions are changed 
accordingly. 

§ 16. Differential Equations of the Second and Higher Orders.— 
Differential equations of the second order may be written in the 


form : : 
dy y 
wa s(s, Y #1) 


Let us introduce the radius of curvature instead of the second 
differential coefficient. Suppose we pass along a curve that 
satisfies the equation and the direction of our motion is deter- 
mined by the angle a it forms with the positive axis of x (counted 
in the usual way from the positive axis of x through ninety 
degrees to the positive axis of y and so on), s being the length of 
the curve counted from a certain point from which we start. 
We then have 


ow tg a, COS @. 
Consequently 
2y Ae da «at da 
dz? cosa dx  cossa ds’ 
or 
a = cos® ast 


da/ds measures the “curvature,” the rate of change of direction 
as we pass along the curve, counted positive when the change 
takes place to the side of greater values of a (if the positive axis 
of x is drawn to the right and the positive axis of y upwards a 
positive value of da/ds means that the path turns to the left). 
Let us count the radius of curvature with the same sign as da/ds 


and let us denote it by p. Then we have 
d 
aioe cos’ a f(a, y, tg a). 


Thus the differential equation of the second order may be said 
to give the radius of curvature as a function of z, y, a, that is to 
say, as a function of place and direction. 


7 
: 
4 

‘ 
5 
- 


——_——s 


Let us assume that this function of three variables is repre- 
sented by a diagram, so that the length and sign of p may quickly 
be obtained for any point and any direction. 

Starting from any given point in any given direction we can 
then approximate the curve satisfying the differential equation 
by a series of circular arcs. Let A (Fig. 92) be the starting point. 
We make M,4 perpendicular to the given direction and equal to 
_ pinilength. For positive values of p, M, must be on the positive 

side of the given direction, for negative 
values on the negative side. M, is 
the center of curvature for the curve 
at A. With M, as center and M,A 
as radius we draw a circular are AB 
and draw the line BM,. On this line 
or on its production we mark the 
point M, at a distance from B equal 
to the value of p that corresponds to 
B and to the direction in which the Fria. 92. 
circular arc reaches B. With M, as 
center and M;,B as radius we draw a circular arc BC and so on. 
The true curve changes its radius of curvature continuously, 
while our approximation changes it abruptly at the points 
A, B,C, ---. The smaller the circular arcs the less will accu- 
rately-drawn circular arcs deviate from the curve. But it must 
be kept in mind that small errors cannot be avoided, when 
passing from one arc to the next. Hence, if the arcs are taken 
very small so that their number for a given length of curve 
increases unduly, the accuracy will not be greater than with 
somewhat longer arcs. The best length cannot well be defined 
mathematically; it must be left to the experience of the draughts- 


yi 4 
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man. 

Some advantage may be gained by letting the centers and the 
radii of the circular arcs deviate from the stated values. The 
circular are AB (Fig. 92) is evidently drawn with too small a 
radius because the radius of the curve increases towards B. If 
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we had taken the radius equal to M,B it would have been too 
large. A better approximation is evidently obtained by making ~ 
the radius of the first circular arc equal to the mean of M,A and 
M,B, and the direction with which it reaches B will also be closer 
to the right direction. 

To facilitate the plotting an instrument may be used consisting — 
of a flat ruler with a hole on one end for a pencil or a capillary - 
tube or any other device for tracing a line. A straight line 
with a scale is marked along the middle of the ruler and a little 
tripod of sewing needles is placed with one foot on the line and 
two feet on the paper. Thus the pencil traces a circular are. 
When the radius is changed, the ruler is held in its position by 
pressing it against the paper until the tripod is moved to a new 
position. By this device the pencil must continue its path in 
exactly the same direction, while with the use of ordinary com- — 
passes it is not easy to avoid a slight break in the curve at the 
joint of two circular arcs. 

Another method consists in a generalization of the method 
for the graphical solution of a differential equation of the first 
order. 

A differential equation of the second order 


CUS ( dy 
de J\ XY; dz 


may be written in the form of two simultaneous equations of the 
first order: 


dy _ 

dx” 

dz 
da 96% Y 2) 


Let us consider the more general form, in which the differential 
coefficients of two functions y, z of x are given as functions of 
eeu ss 


dy 
de 1 Y 2); 


=e 
4 
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dz 
dx ma) (x, Y, 2). 


We may interpret 2, y, z as the codrdinates of a point in space — 
and the differential equation as a law establishing a certain 
direction or the opposite at every point in space where f(z, y, 2) 
and g(x, y, 2) are defined. A curve in space satisfies the dif- 
ferential equation, when it never deviates from the prescribed 
direction. Its projection in the zy plane represents the function 
y and its projection in the zz plane represents the function z. 

Let us represent y and z as ordinates and 2 as abscissa in the 
same plane with the same system of codrdinates. Any point in 


Fie. 93. 


space is represented by two points with the same abscissa. The 
functions f(x, y, 2) and g(z, y, z) we suppose to be given either 
by diagrams or by certain methods of construction or calculation. 
For any point that we have to deal with, the values of f(z, y, 2) 
and g(z, y, 2) are plotted as ordinates to the abscissa x, but for 
clearness sake not in the same system of codrdinates as y and 2, 
but in another system with the same axis of ordinates and an 
axis of x parallel to the first and removed far enough so that the 
drawings in the two systems do not interfere with one another. 
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Starting from a certain point P(2p, yp, Zp) in space we represent 
it by the two points Pi(ap, yp) and P2(ap, z,) in the first system 
and the values of f(xp, yp, Zp) and g(xp, Yp, Zp) by the two points 
A, and Az in the second system of coérdinates (Fig. 93). The 
points A; and A, determine certain directions MAy;’, and MA,’ 
of the curves 2, y and 2, z, the point M (Fig. 93) being placed at a 
distance from the axis of ordinates equal to the unit of length by 
which the ordinates representing f(z, y, z) and g(a, y, 2) are 


measured. Through P; and P, we draw parallels to MA,’ and © 


MA,’ as far as Q; and Q. with the codrdinates aq, yq and 2g, 2. 


With these coérdinates the values f(a, Yq; %q) and g(%q; Yay 2a) 
are determined, which we represent by the ordinates of the © 


points B,, B,. These points again determine certain directions 
‘ parallel to which the lines Qi; and Q.R, are drawn, etc. In this 
manner we find first approximations y; and % for the functions 
y and z and corresponding to these approximations we find 
curves representing f(x, y1, 21) and g(x, y1, 21). These curves are 
now integrated graphically, the integral curve of f(z, ym, %) 
beginning at P, and the integral curve of g(x, m1, 2) at Pe and 
lead to second approximations y and z: 


Yo = Yo a it Ac? Ys z)dz, 


2 = Sp © jp g(x, 1; z)dax. 


For these second approximations the values of f(a, ye, 2) and 
g(&, Y2, 22) are determined at a number of points along the curves 
2, y and x, % sufficiently close to construct the curves representing 
F(x, Ye, %) and g(x, y, %). By their integration a third approxi- 
mation y3, 23 is obtained 


Ys = Yo + { f(x, Yo, x) da, 


3 = Sp zie ) g(x, Y, z)dx, 
*p 


j a 
~ 
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and so on as long as a deviation of an approximation from the 
one before can still be detected. As soon as there is no deviation 
for a certain distance z — x, the curve represents the true solu- 
tion (as far as the accuracy of the drawing goes). The curve is 
continued by taking its last point as a new starting point for a 
similar operation. 

The distance over which the intégral is taken can in general 
not surpass a certain limit where the convergence of the approxi- 
mations ceases. But we are free to make it as small as we please 
and accordingly increase the number of operations to reach a 
given distance. It is evidently not economical to maké it too 
small. On the contrary, we shall choose it as large as possible 
without unduly increasing the number of approximations. 

In the case of a differential equation 


a a 
ae, 2 Ys ae 


we have f(a, y, 2) = z, and the curve z, 2 is identical with the 
curve representing the values of f(z, y, z). We shall therefore 
draw it only once. 

The proof of the convergence of the approximations is almost 
the same as in the case of the differential equation of the first 
order. 

For the n + 1* approximation we have 


Yntt = Yor HEE (2, Yn, Zn)dt; 2np1 = Zp + ih ; G(X, Yny Zn) dx. 
For the true ike that passes through the Sane Ln, Ups Sp We 
- find by integration 

— +f fe y, 2)dx; 2= at [ole Y, 2)des 
hence y : 


Yui Y = f ue, Ynys Zn) Es Ste, UP 2) |dx; 


P 


amis — 2 =f [ales Yas to) — 90 vy, 2) 
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Now let us write 


ya, Ynys Zn) _ gi ecp Y, 2) = F(&, Yny Bn) — fs Ys Zn) 


(Yn — y) 


Yn—Y 
X,Y, Zn) — J\X, Y;, % 
flint) = feng, — 9 
and similarly 
Ly Ynys on) — L,Y; Sn 
g(x, Yn Zn) — g(a, Y; z) = Ges Yn) ee ee y) 
Yn — Y 
ae g(x, Yy; Zn) a g(x, Y; z) (Zn _ z). 


Zn — 
The quotients of differences 
f(z, Yns Zn) iam f(a, Ys Zn) 
Yn — Y 

and the three others are equal to certain values of df/dy, df/d2, 
dg/dy, Og/dz for values of y, 2 between y and y, and between 2 
and Zn (Y, Yn, 2, Zn not excluded). Let us assume that for the 
region of all the values of 2, y, z concerned the absolute value of 
of/dy and of/dz, is not greater than M,, and that of dg/dy and 
dg/dz not greater than Ms, and that 6, e, denote the maximum 
of the absolute values of y — y, and z— 2, in the interval 
a, toa. Then it follows that the absolute values of 


ite Yny Zn) ee F(a, Y, 2) and g(x, Yns Zn) pe g(x, Y; 2) 
are not greater than 


Mi(6n — En) and Mo(6n + En). 


Hence for the maximum values of yny1 — y and Zn41 — 2, which 
are denoted by 6,41 and €n4; we obtain the limits 


On4t = Mi(6,, Sis En) | X— Lp | ’ Ent1 < Mo(6n oe En) | t— Lp | ’ 


and 
Ont =e Enti = (M, == Mz) | vL— Lp | (6n + ry 


If therefore the interval x — x, of the integration is so far 
reduced that 


(Mi + Me)|a —a, Sk <1, 


+= Mi 


el i i 
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dntt + é€nti is not larger than the fraction k of (8, + en), but 


_ from the same reason 


(Sn ‘17 En) S K(Sn—1 st Eni) (6n1 2 €n—1) = k(Sn-2 ae En—2), etc. ; 


therefore 
Bn4t + éengr S h"(S + a). 


That is to say, for a sufficiently large value of m 5n41 and eny1 
will both become as small as we please. 

As in the case of the differential equation of the first order it is 
not worth while, as a rule, to investigate the convergence for the 
purpose of finding a sufficiently close approximation by graphical 
methods. It is-better at once to tackle the task of drawing the 
approximations and to repeat the operations until no further 
improvement is obtained. The curve will then satisfy the 
differential equation as far as the graphical methods allow it 
to be recognized. 

When the values of f(z, y, z) or g(x, y, 2) become too large 
we can have recourse to the same device that we found useful 
with the differential equation of the first order. Instead of 2, 
one of the other two variables y or z may be considered as inde- 
pendent, so that the equations take the form 


de _iil dz _ gl, y, 2) 

dy f(z, y,2)’ dy f(x,y,2)’ 
or 

dx 1 dy _ f(x,y, 2) 


dz “i g(a, y; Bye dz - g(x, Y; ae 


or we may introduce a new system of codrdinates 2’, y’, 2’ and 
consider the resulting differential equations. 

The second method for the integration of differential equations 
of the first order can also be generalized to include the second 
order. Let us again consider the more general case 


dy. dz 
a = f(a, Y; z), dz ad g(x, Y; Z). 


Starting from a point 2, y, 2 the changes of y and z (denoted by 
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; a 
k and 1) can be calculated for a small change h of x by the fol-— 
lowing formulas analogous to those used for one differential — 


equation of the first order: 


hy = f(z, Y, z)h; = g(x, Y; zh; 


h k l h k | 
ke =i (+5, ot oe +8 hs b=g(2+5, yt, 2+ yh | 


lL. h k 
ks =f(et}, if z+ 3) b= (ety, Te +5 Ji 
kg=f(@ th, yth, z+h)h; L=gath yt+kh, 2+ bh; 


_ kth ha ke | ,o bth ,_ be 
a 2 re! a 9 2 fae 2 Pa! ip imme 2 5) 


and with a high degree of approximation, 
b= pr 3g — fp); b= p a 


These calculations may be performed graphically. For this 
purpose the functions f(x, y, z) and g(a, y, 2) must be given in 
some handy form. We notice that in our formulas the first 
argument assumes the values 2, x + h/2, x+h. In the next 
step where x + h, y + k,2-+ lare the coérdinates of the starting 
point that play the same part that 2, y, z played in the first 
step, we are free to make the change of the first argument the 
same as in the first step, so that in the formulas of the second 
step it assumes the values x + h, e+ $h, x + 2h and so on for 
the following steps. All the values of the first argument can 
thus be assumed equidistant. Let us denote these equidistant 
values by 


U0, V1, V2, Uz, ° °° 


The values of f(, y, 2) and g(2, y, z) appear in all our formulas 
only for the constant values 


v= Xo, M1, Xa, °° 


For each of these constants f and g are functions of two inde- 
pendent variables and as such may be represented graphically 


hs ie 
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by drawings giving the curves f = const. and g = const., each 
value of x corresponding to a separate drawing. These drawings 
we must consider as the graphical form in which the differential 
equations are given. It may of course sometimes be very tire- 
some to translate the analytical form of a differential equation 
into a graphical form, but this trouble ought not to be laid to 
the account of the graphical method. 

The method now is similar to that used for the differential 
equation of the first order. y and z are plotted as ordinates in 
the same system in which z is the abscissa. Equidistant parallels 
to the axis of ordinates are drawn 


x= X%, © = M1, T= M, ete. 


- On the first x = x we mark two points with ordinates yo and 2p, 


and from the drawing that gives the values of f(xo, y, 2) and 


iS : 
<--> 


— 


' 
1 


Fig. 94. 


g(xo, y, 2) as functions of y and z we read the values f(xo, yo, 20) 
and g(2o, Yo, 20) and draw the lines from ao, yo, and 2, % to the 
points 

tm, Yoth, and m, 2+h. 
The intersections of these lines with the parallel 2 = 2, furnishes 


the points 
11 
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/ l 
MH, yo+ "5 and 1, 20+ 5. 


With these ordinates we find from the second drawing the values 
k l k l, 
f (+3, at3) and a( 2 w+ Ayoth), 


and by their help we can draw the lines from 2, Yo and 2, 20 


to the points 
X25 Yo + ko and a, yo+ bh. 


The intersections of these lines with the line 2 = 2; furnishes the 
points 


k l 
ets yot > and 7, Zo + oe 


and with these ordinates we find the values 


k. ly ko /} 
f(a mt 2043), (a, w+S,0+2), 


which enable us to draw the lines from 2, yo and 2a, 2 to 2, 
yo + kz and a, 20 + 3. 

With these two ordinates we find from the third diagram (2 = a) 
the values 


f (2, Yo + ks, 20+ 13) and g (x2, Yo + kg, 20 + ls), 


which finally enable us to draw the lines from zy and 292% to 
X2, Yo + kg and x, 29 + ly. 

On the vertical line x = a, we thus obtain four points, By, Bs, 
Bz, Bs, corresponding to yo + ki, yo + he, yo + ks, yo + ka and 
four points, By’, B,’, B;’, Ba’, corresponding to zo + h, zo + bh, 
zo + ls, 20 + ly (Fig. 94). | 

B.B3 and B,B, are bisected by the points C; and C; By’ B,’ 
and By’B,' by the points Cy’, C2’. Finally C,C, and C1’/C,’ are 
divided into three equal parts and the points B and B’ are found 
in the dividing points nearest to C; and Cy’. 

The same construction is then repeated with B and B’ as 
starting points and furnishes two new points on the vertical 


———— 


ie 
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x = x, and so on. To test the accuracy the construction is 
repeated with intervals of x of double the size. The difference 
in the values of y and of z found for 2 = 2 enables us to estimate 
the errors of the first construction—they are about one-fifteenth 
of the observed differences. 

Both methods are without difficulty generalized for the integra- 
tion of differential equations of any order. We can write a 
differential equation of the nth order in the form 


d"x dx ie 
ass aah ee =) 


or in the form of n simultaneous equations of the first order 


dx 
dt = %; 
dx, 
Pip 
dxtn-2 
di Un-1; 
Gent 


dt = 7 Uv, U1; Yo, °° * Ln-1). 


A more general and more symmetrical form is 
dx 
di = AG ©, U1, °° * Ln—1)s 


dx 
oe = falt, @,%;,°°° ina) 


din- 
dt : = Fall, XV, U1, °°° Ent): 


The functions 2, 21, #2, **+ n-1 are then represented as ordinates 
to the abscissa t, so that we have n different curves. When the 
function f(t, 2, 21, %, °** %n-1) is given ina handy form, so that 
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its value may be quickly found for any given values of ¢, x, a, 
.++ a@,-1, there is no difficulty in constructing n curves whose 
ordinates represent the functions 2, a1, %, **: %4. Starting 
from given values of t, 2, 21, %2, +++ &n-1 we have only to apply 
the same methods that have been explained for the first and the 
second order. 
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